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Abstract 

We have studied hadronic events from e + e~ annihilation data at centre-of-mass energies 
of y/s =172, 183 and 189 GeV. The total integrated luminosity of the three samples, 
measured with the OPAL detector, corresponds to 250 pb _1 . 

We present distributions of event shape variables, charged particle multiplicity and 
momentum, measured separately in the three data samples. From these we extract mea- 
surements of the strong coupling a s , the mean charged particle multiplicity (n ch ) and the 
peak position £ in the £ p = ln(l/x p ) distribution. 

In general the data are described well by analytic QCD calculations and Monte Carlo 
models. Our measured values of a s , (n c h) and £o are consistent with previous determina- 
tions at \/s=M z o. 
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1 Introduction 



We study the general features of hadronic decays in e + e~ — > (Z /^)* — > qq reactions 
at the highest available centre-of-mass (cm.) energies, as a continuation of our earlier 
publications at cm. energies of y/s = 130 — 136 GeV [I] and y/s = 161 GeV 0. 

In the last three years LEP produced e + e~ collisions at centre-of-mass energies of 
y/s ~ 172, 183 and 189 GeV. The total integrated luminosity measured with the OPAL 
detector at these energies corresponds to approximately 250 pb _1 . 

Previous studies using e + e~ annihilation data at cm. energies up to y/s = 183 GeV 
have shown that QCD based models and calculations give a good description of the 
observed data fl]-|10|. Here we examine the overall consistency of QCD at yet higher cm. 
energies, with improved statistical precision owing to the dramatic increase in integrated 
luminosity at cm. energies above y/s = 161 GeV. 

We determine the strong coupling strength a s at each of the three energies. The 
large data sample at 189 GeV, together with the fact that the hadronization corrections 
become smaller at higher cm. energies, allows a precision comparable with that achieved 
at y/s = M z o using event shapes. We compare our results to those obtained at lower cm. 
energies, notably at y/s = M z o. We also determine the charged particle distributions and 
study the evolution with the cm. energy of the mean charged particle multiplicity and 
the peak position £ in the £ p = ln(l/x p ) distribution. Jet-multiplicity related observables 
have been analized and used to determine a s at these cm. energies in a separate study 
of e + e~ annihilation data taken by the JADE and OPAL experiments at cm. energies 



between 35 and 189 GeV 11 



The majority of hadronic events produced at cm. energies above the Z° resonance 
are 'radiative' events in which initial state photon radiation reduces the energy of the 
hadronic system to about M z o. An experimental separation between radiative and non- 
radiative events is therefore required. In addition, at energies above 160 GeV, W + W~ 
production becomes kinematically possible and contributes a significant background to 
the (Z°/7)* — > qq events. For y/s > 180 GeV, Z°Z° production must also be taken into 
account. 

In this paper we use similar techniques to those of our previous analysis of e + e~ 
annihilation data at cm. energies of 130-136 GeV [0 and 161 GeV ||. In Section ||] we 
briefly describe the OPAL detector, in Section [3] we describe the samples of data and 
simulated events and in Section [| we explain our selection requirements and analysis 
procedures. In Section |5| we then present the results in the form of event shape variables, 
determination of a s , charged particle multiplicities and momentum spectra. Section |] 
gives a summary and conclusion. 



2 The OPAL detector 

The OPAL detector operates at the LEP e + e~ collider at CERN. A detailed description 



can be found in reference ||12|| . The analysis presented here relies mainly on the mea- 
surements of momenta and directions of charged tracks in the tracking chambers and of 
energy deposited in the electromagnetic and hadronic calorimeters of the detector. 
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All tracking systems are located inside a solenoidal magnet which provides a uniform 
axial magnetic field of 0.435 T along the beam axisQ. The magnet is surrounded by a lead 
glass electromagnetic calorimeter and a hadron calorimeter of the sampling type. Outside 
the hadron calorimeter, the detector is surrounded by a system of muon chambers. There 
are similar layers of detectors in the forward and backward endcaps. 

The main tracking detector is the central jet chamber. This device is approximately 
4 m long and has an outer radius of about 1.85 m. It has 24 sectors with radial planes of 
159 sense wires spaced by 1 cm. The momenta p of tracks in the x-y plane are measured 
with a precision a p /p = yll.02 2 + (0.0015 ■ p[GeV/c]) 2 . 

The electromagnetic calorimeters in the barrel and the endcap sections of the detector 
consist of 11704 lead glass blocks with a depth of 24.6 radiation lengths in the barrel and 
more than 22 radiation lengths in the endcaps. 



3 Data and Monte Carlo samples 

The three data samples at ^/s ~ 172, 183 and 189 GeV that are used in this analysis 
were recorded as part of the LEP-2 programme between 1996 and 1998. The luminosities, 
evaluated using small angle Bhabha collisions, and the mean cm. energies are tabulated 
in Table |I[ 

Monte Carlo event samples were generated at cm. energies of 172.0, 183.0 and 189.0 GeV, 
including full simulation of the OPAL detector f!3]|. Events for the process e + e" 



[ZP/'j)* — > qq, referred to as u (Z Q /j)* events" , were generated using the PYTHIA 5.722 |TJ 



parton shower Monte Carlo code with initial and final state photon radiation and frag- 
mentation of the parton final state handled by the routines of JETSET 7.408 [14]. The 
simulation parameters have been tuned to OPAL data taken at the Z° peak ||15|| . The 
JETSET Monte Carlo is able to provide a good description of experimental data for e + e~ 
annihilations with cm. energies from 10 GeV to 161 GeV |l|,|3|-|6|,[T0[|. 

As an alternative fragmentation model, we generated events with the HERWIG 5.9 [0 
parton shower Monte Carlo program, also tuned to OPAL data, as described in H,|T5j|. 



As an alternative model for initial state radiation (ISR) the YFS3ff 3.6 [I7| generator was 
used, coupled to the same JETSET routines to handle the parton final states. 

In addition we generated events of the type e + e~ — > 4 fermions (diagrams without 
intermediate gluons). These 4-fermion events, in particular those with four quarks in 
the final state, constitute a major background for this analysis. Simulated 4-fermion 
events with quarks and leptonic final states were generated using the GRC4F 1.2 Monte 



Carlo Model ||18|| . The final states were produced via s-channel or t-channel diagrams 
and include W + W~ events. This generator is interfaced to JETSET 7.4 using the same 
parameter-set for the parton shower, fragmentation and decays as for (Z°/j)* events. 
Two other possible sources of background events were simulated. Hadronic two-photon 



processes were evaluated using PYTHIA, HERWIG and PHOJET ||19[| . Production of 
e+e" -> (Z°/7)* -> rf was evaluated using KORALZ |EB. 



1 In the OPAL coordinate system the x axis points towards the centre of the LEP ring, the y axis 
points upwards and the z axis points in the direction of the electron beam. The polar angle 9 and the 
azimuthal angle <p are defined w.r.t. z and x, respectively, while r is the distance from the z-axis. 
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In addition to the Monte Carlo event generators discussed above, we use the event 

ARIADNE is used to provide a 



generators ARIADNE 4.08 [0] and COJETS 6.23 |2 



systematic check on the hadronization corrections for our a s measurement (Section |5.2|) . 
In contrast to the the other models, COJETS is based on a parton shower model using 
independent fragmentation and does not take coherence effects into account. Both CO- 
JETS and ARIADNE are used in addition to PYTHIA and HERWIG to compare to our 
corrected data distributions. The parameter sets used for ARIADNE and COJETS are 
documented in [[H]] and [p3|] ; both models provide a good description of global e + e~ event 
properties at y/s = M z o, as do PYTHIA and HERWIG. 



4 Analysis method 
4.1 Selection of events 

Hadronic events are identified using criteria as described in ||24|| . The efficiency of se- 
lecting non-radiative hadronic events is essentially unchanged with respect to lower cm. 
energies and is approximately 98% [p5fl . We define as particles tracks recorded in the 
tracking chambers and clusters recorded in the electromagnetic calorimeter. The tracks 
are required to have transverse momentum to the beam axis px > 150 MeV/c, a number 
of hits in the jet chamber Nhus > 40, a distance of the point of closest approach to the 
collision point in the r — plane do < 2 cm and in the z direction zq < 25 cm. The clusters 
in the electromagnetic calorimeter are required to have a minimum energy of 100 MeV in 
the barrel and 250 MeV in the endcap sections J^j. 

To reject background from e + e~ — > t + t~ and 77 — > qq events and to ensure the events 
are well contained in the OPAL detector we require at least seven accepted tracks, and 
the cosine of the polar angle of the thrust axis | cos 6t\ < 0.9. The number of events after 
this preselection, for each cm. energy, are listed in Table [1]. 

To reject radiative events, we determine the effective cm. energy y/s' of the observed 
hadronic system as follows pEfl. Isolated photons in the electromagnetic calorimeter are 



identified, and the remaining particles are formed into jets using the Durham |27j algo- 
rithm with a value for the resolution parameter y cnt = 0.02. The energy of additional 
photons emitted close to the beam directions is estimated by performing three separate 
kinematic fits assuming zero, one or two such photons. The probabilities of the fits are 
used to select the most likely of these three. The value of y/s* is then computed from 
the fitted momenta of the jets, excluding photons identified in the detector or close to 
the beam directions. The value of \fs' is set to \/s if the fit assuming zero initial state 
photons was selected. The 4-momenta of the measured particles are then boosted into 
the rest frame of the observed hadronic system. 

To reject events with large initial-state radiation (ISR), we require \/s — yfs' < 10 GeV. 
This is referred to as the 'ISR- fit' selection. In Figure |I|a we compare the xTs* distribution 
of the data to simulated (Z /^)* and 4-fermion events. The simulated (Z°/7)* events are 
classified into radiative events with y/s — y/s' tTUC > 1 GeV (where \/ s' trne is the true effective 
cm. energy, determined from Monte Carlo information), and non-radiative events, which 
is the complement. While about 27% of the selected (Z°/ / y)* events are by this definition 
radiative events, the fraction of events with yfs — \J s' tvuc > 10 GeV is only 5%. The 
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contributions of 4-fermion events, simulated with GRC4F, are indicated separately. The 
background from 4-fermion events and the efficiency of selecting non-radiative events are 
given in Table [I]. 

The background from e + e~ — > t + t~ and two-photon events of the type 77 — > qq is 
estimated from Monte Carlo samples to be less than 0.3% and is neglected. 

In order to reduce the background of 4-fermion events on the remaining sample, we 
test the compatibility of the events with QCD-like production processes. A QCD event 
weight Wqcd is computed as follows. We force each event into a four-jet configuration in 
the Durham jet scheme and use the EVENT2 [28 program to calculate the 0{a 2 s ) matrix 
element \M(pi, p 2 , P3, p^f for the processes e + e~ — > qqqq, qqgg pSfl . Since neither quark 
nor gluon identification is performed on the jets, we calculate the matrix element for each 
permutation of the jet-momenta and use the permutation with the largest value for the 
matrix element to define the event weight: 

Wqcd = max log (\M(pi,p 2 ,P3,Pa)\ 2 ) (1) 

with pi the momenta of reconstructed jets. Note that the definition of the event weight 
Wqcd contains kinematic information only and is independent of the value of a s . 

This weight Wqcd is expected to have large values for processes described by the QCD 
matrix element, originating from (Z /^)* — > qq, and smaller values for W + W~ events. In 
Figure [L]b we compare the data distribution of Wqcd to the expectations of our simulation. 
A clear separation between the (Z /^)* and 4-fermion events is achieved by requiring 
Wqcd > —0.5. This requirement reduces the 4-fermion background considerably, as can 
be seen in Table [l], whereas the efficiency for selecting non-radiative (Z°/7)* events is 
only slightly reduced. In Figure we show, as an example, the effect of the 4-fermion 



background rejection on the distribution of thrust T, as defined in Section |5.1| . The 
background from hadronic decays of W + W~ events has predominantly low values for T, 
and is largely removed by our selection. 



4.2 Correction procedure 

The remaining 4-fermion background in each bin has been estimated by Monte Carlo 
simulation and is subtracted from the observed bin content. A bin-by-bin multiplication 
procedure is used to correct the observed distributions for the effects of detector resolution 
and acceptance as well as for the presence of remaining radiative (Z /^)* events. For the 
multiplicity measurement presented in Section |5.3|, we use a matrix correction procedure 



to account for detector effects rather than a bin-by-bin procedure. 

For the bin-by-bin correction procedure, each bin of each observable is corrected from 
the "detector level" to the "hadron level" using two samples of Monte Carlo (Z /^)* 
events at each cm. energy. The hadron level does not include initial state radiation or 
detector effects and allows all particles with lifetimes shorter than 3 x 10 -10 s to decay. 
The detector level includes full simulation of the OPAL detector and initial state radiation 
and contains only those events which pass the same cuts as are applied to the data. The 
bin-by-bin correction factors are derived from the ratio of the distributions at the hadron 
level to those at the detector level. 
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A bin-by-bin correction procedure is suitable for most quantities as the effects of finite 
resolution and acceptance do not cause significant migration (and therefore correlation) 
between bins. For the multiplicity measurement however, such a method is not readily 
applicable, due to the large correlations between bins. 

The four-momenta of tracks and of the electromagnetic calorimeter clusters not asso- 
ciated with tracks were used to calculate event shapes. When a calorimeter cluster had 
associated tracks, their expected energy deposition was used to reduce double counting 
by correcting the cluster energy. If the energy of a cluster was smaller than the expected 
energy response of the associated tracks, the cluster energy was not used. The masses 
of all charged particles were set to the pion mass and the invariant masses of the energy 
clusters were assumed to be zero. 

4.3 Systematic uncertainties 

The experimental systematic uncertainty is estimated by repeating the analysis with var- 
ied experimental conditions. In order to reduce bin-to-bin fluctuations in the magnitudes 
of the systematic uncertainties we average the relative uncertainty over three neighbouring 
bins. 

To allow for any inconsitencies caused by differences between the responses of the 
tracking or the calorimeter, three differences are formed for the quantities in each bin; 
between the standard result, the one obtained using tracks and all clusters and the one 
obtained using only tracks. The largest of them is taken as the systematic error. 

The inhomogeneity of the response of the detector in the endcap region was allowed 
for by restricting the analysis to the barrel region of the detector, requiring the thrust 
axis of accepted events to lie within the range | cos#t| < 0.7. The event sample is reduced 
in size by approximately 27%. The corresponding systematic error is the deviation of the 
results from those of the standard analysis. 

For observables measured using information from charged particles only, we evaluate 
an additional uncertainty on the track modelling. The maximum allowed distance of the 
point of closest approach of a track to the collision point in the r-<p plane do was changed 
from 2 to 5 cm, the maximal distance in the z direction zq from 25 to 10 cm and the 
minimal number of hits from 40 to 80. These modifications change the number of tracks 
by up to approximately 12%. The quadratic sum of the deviations from the standard 
result is taken to be the systematic error due to the uncertainty of the track modelling. 

Uncertainties arising from the selection of non-radiative events are estimated by re- 
peating the analysis using a different technique to determine the value for y/s' \ This 
technique differs from our standard y/s' algorithm in that in this case the kinematic fit 
assumes always one unobserved photon close to the beam direction for each event. The 
final event sample with this y/s' algorithm has an overlap of approximately 94% with the 
standard sample, and is reduced in size by 3%. The difference relative to the standard 
result is taken as the systematic error. 

We evaluated our (Z°/j)* selection efficiencies using events generated with YFS3ff, 
which contains QED exponentiated matrix elements for ISR up to 0(a% M ). Similar esti- 
mates of efficiencies are obtained if this model is used instead of PYTHIA, and therefore 
no additional systematic uncertainty is assigned. 
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Systematic uncertainties associated with the subtraction of the 4-fermion background 
events are estimated by varying the position of the cut-value on the QCD event weight 
Wqcd (see Figure |I|b). We vary this position to Wqcd > —0.8 which increases the 
event sample by approximately 7%, and Wqcd > which decreases the event sample by 
approximately 12%. We take the maximum deviation from our standard result Wqcb > 
—0.5 as the systematic uncertainty. 

In addition we vary the predicted background to be subtracted, within its measured 
uncertainty of 5% and use the largest difference from the standard result as a sys- 
tematic uncertainty. 

The difference in the results when we use simulated (Z°/j)* events generated using 
HERWIG instead of PYTHIA is taken as the uncertainty in the modelling of the (Z /^)* 
events. 



5 Results 

5.1 Event shapes 

The properties of hadronic events may be characterised by a set of event shape observables. 
The following quantities are considered: 



Thrust T: defined by the expression pl|,p2 

T = max(%^) . (2) 

The thrust axis ht is the direction n which maximises the expression in parenthesis. 
A plane through the origin and perpendicular to divides the event into two 
hemispheres Hi and H 2 . 

Thrust major T ma j or . The maximisation in equation [2] is performed with the condition 
that n must lie in the plane perpendicular to Hr- The resulting vector is called 
fir ■ ■ 

major 

Thrust minor T m i nor . The expression in parenthesis is evaluated for the vector ^T minor 
which is perpendicular both to % and to % . . 

IT IT J- J maior 



Oblateness O. This observable is defined by O = T ma j or — T minor p3 | 



Sphericity S and Aplanarity A. These observables are based on the momentum ten- 
sor 



b H = — — — , a,(J=l,2,3 



The three eigenvalues Qj of S a/3 are ordered such that Qi < Q 2 < Qz- These then 
define S pl,03] and A p6[ by 



S = ^(Qi + Q 2 ) and A=^Qi 
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C-parameter. The momentum tensor S a " is linearised to become 

ao/J YU(P?Pi)/\Pi\ „ a i o Q 

B p = — ; ; , a,/* =1,2,3 . 

J2i \Pi\ 

The three eigenvalues Aj of this tensor define C [37] with 



C = 3(AiA 2 + A 2 A3 + A 3 A 1 ) . 
Heavy Jet Mass Mh- The hemisphere invariant masses are calculated using the parti- 



cles in the two hemispheres H 1 and H 2 . We define Mh |p8| , [39|] as the heavier mass, 
divided by yjs . 

Jet Broadening variables B T and B\y. These are defined by computing the quantity 

[EieH k \Pi x n T \\ 



Bi 



2E, \Pi 



for each of the two event hemispheres, Hk, defined above. The two observables [|4(J 
are defined by 

Bt — B\ + B2 and Byy = max(5i, B2) 
where Bt is the total and By/ is the wide jet broadening. 

Transition value between 2 and 3 jets y® 3 . The value of y cnt for the Durham jet 
scheme at which the event makes a transition between a 2-jet and a 3-jet assignment. 

In the following, we use the symbol y to denote a generic event shape observable, 
where larger values of y indicate regions dominated by the radiation of hard gluons and 
small values of y indicate the region influenced by multiple soft gluon radiation. Note that 
thrust T forms an exception to this rule, as the value of T reaches one for events consisting 
of two collimated, back-to-back, jets. We therefore occasionally use 1 — T instead. 

Figures || and [| show the distributions of the event shape observables T, T ma j or , T minor , 
A, C, M H , S, O, Bt, B w , and y^ 3 , corrected for detector acceptance and initial state 
radiation, plotted at the weighted centres of the bins. The data obtained at 189 GeV 
are shown with the statistical and systematic uncertainties added in quadrature. The 
numerical values for all the event shape distributions, for cm. energies of 172, 183 and 
189 GeV, are listed in Tables g to g. 

With the statistics available at 189 GeV, all event shapes are determined with high 
precision. There are no deviations from the predictions of PYTHIA, ARIADNE and 
HERWIG. The COJETS Monte Carlo model deviates from the data for T minor , Mh and 
O, but describes the remaining event shapes reasonably well. Note that all generators 
have been tuned at ^/s = M z o, and are not re-tuned at these higher energies. 

The mean value of the thrust distribution, (T), as a function of the cm. energy, is 
shown in Figure [| together with data from lower energy measurements [|l],0,|25|,|41| and the 



predictions for the energy evolution from PYTHIA, HERWIG, ARIADNE and COJETS. 
The predictions of all four Monte Carlo models are consistent with our new measurements. 
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5.2 Determination of a. 



s 



Our measurement of the strong coupling strength a s (Q) is based on fits of the QCD 
predictions to the corrected distributions for 1 — T, M H , C, B T , B w and y^- The 
theoretical descriptions of these six observables are the most complete, allowing the use of 



combined 0(a^)+NLLA QCD calculations ||40| , p5| -|45H . We follow the procedures described 



in references |],^| as closely as possible in order to obtain results which we can compare 
directly to our previous analysis. In particular, we choose the so-called In (R) -matching 
scheme, and fix the renormalization scale parameter, = fi/y/s, to x M = 1, where /i is 
the energy scale at which the theory is renormalised. The (9(a^)+NLLA prediction for C 
has only recently become available ||45|| , and was not used in our earlier publications. 



Analytic QCD predictions describe distributions at the level of quarks and gluons 
(parton level). The predictions are convolved with hadronization effects by multiplying 
by the ratio of hadron- to parton level distributions determined using a Monte Carlo 
model. We use PYTHIA to generate events at yfs = 172, 183 and 189 GeV for this 
purpose. 

The fit ranges are determined by the following considerations. For each observable, the 
ratio between Monte Carlo distributions computed for partons and hadrons is required to 
be unity to within about 10% and the distribution of partons from Monte Carlo models 
is required to be described well by the analytic predictions. The fit ranges are identical 
to the fit ranges used in our studies at ^/s = M z o, 130 and 161 GeV fl|,p|,[2~5|l. 

We find satisfactory fits for all six observables, for all three cm. energies. In particular, 
the C distributions are well described by the newly introduced predictions. In Table [5] we 
quote the result on a s (172 GeV), a s (183 GeV) and a s (189 GeV) for all six observables. 
In Figure |] we show the results for the determination of a s (189 GeV). Note also that 
outside the fit ranges the data agree with the 0(a;?)-|-NLLA predictions. 

For each observable, the statistical uncertainties are estimated from the variance of a s 
values derived from fits to 100 independent sets of simulated events, each with the same 
number of events as the data. 

We also derive a combined result of the six observables for the strong coupling strength 



using the weighted average procedure as described in reference ||25|| , which includes the 
correlations between the shapes. The statistical uncertainty of the combined result is 
estimated using 100 independent samples of Monte Carlo events in the same manner as 
for the individual measurements. 

The experimental uncertainty is estimated by adding in quadrature the following con- 
tributions: the largest difference between the central result and the results from tracks 
and all clusters or tracks only; the difference found when using the alternative VV selec- 
tion; the difference when requiring the thrust axis to lie in the range | cos^t] < 0.7; the 
difference from using the variation on IVqcd to reject 4-fermion events; and the difference 
when the 4-fermion background is scaled by ±5%. The experimental uncertainties, shown 
for each observable in Table |5], are of similar size as the statistical uncertainties. In Table |] 
we present details of the experimental uncertainty for the weighted mean value of a s , at 
each cm. energy. 

The hadronization uncertainty is defined by adding in quadrature: the larger of the 
changes in a s observed when varying the hadronization parameters b and ctq by ± 1 
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standard deviation about their tuned values [15| in PYTHIA; the change observed when 



the parton virtuality cutoff parameter is altered from Q = 1.9 GeV to Q = 4 GeV 
in PYTHIA (without changing the other Monte Carlo parameters), corresponding to a 
change of the mean parton multiplicity from 7.4 to 5.0; the change observed when at 
the parton level only the light quarks u, d, s and c are considered in order to estimate 
potential quark mass effects; and both differences with respect to the standard result 
when HERWIG or ARIADNE are used to account for the hadronization effects, rather 
than PYTHIA. For all observables the hadronization uncertainties are given in Table |^. 
These uncertainties are relatively small compared to the statistical uncertainty of a s for 
the 172 GeV sample, but become of similar size to the statistical uncertainty for the 
189 GeV sample. In Table ||| we show the individual contributions to the hadronization 
uncertainty for the weighted mean value of a s . 

The importance of uncomputed higher order terms in the theory may be estimated 
by studying the effects of varying the renormalization scale parameter x M . We estimate 
the dependence of our fit results on the renormalization scale x M as in references |l,|25 



by repeating the fits using x^ = 0.5 and = 2. We define the average deviation from 
the central result as the systematic uncertainty. We find variations which are generally 
larger than any other systematic variation and which are highly correlated between all 
observable^. 

The total uncertainty for each individual observable is computed by adding in quadra- 
ture the statistical, experimental, hadronization and scale uncertainties. The total uncer- 
tainty on a s is typically around 10% at y/s = 172 GeV and around 5% at \/s = 183 and 
y/s = 189 GeV. 

As a consistency check we repeated the fits to a s while varying the fit ranges. The 
fit range was changed by plus or minus one bin at either side. In general the observed 
deviations to the central value of a s were small and we do not include this check in our 
systematic uncertainty. 

Our result for the computed average of the six event shapes, each weighted with its 
total uncertainty, is 

a s (172 GeV) = 0.092 ± 0.006(stat.) ± 0.008(syst.), 

a s (183 GeV) = 0.106 ± 0.003(stat.) ± 0.004(syst.), (3) 

a s (189 GeV) = 0.107 ± 0.001(stat.) ± 0.004(syst.). 

The systematic uncertainty has contributions from experimental effects ±(0.002 — 0.007), 
hadronization effects ±(0.001 — 0.002) and from variations of the renormalization scale 
±(0.002 — 0.004), as explained above. The systematic variations of the combined results 
are detailed in Table ^| Their values evolved to M^o are listed in Table [7]. In Figure [F] we 
show the values of a s we obtain from all event shapes at the three energies. In the same 
figure we also show the computed weighted averages of a s . The spread in the values of 
a s as obtained from the six event shapes is small, and the computed weighted average 
of a s covers the individual measurements within the uncertainty. Note however that a s 



2 The compensation in a s due to a change of the renormalisation scale in the QCD predictions is 
proportional to the value of a s itself. The scale uncertainty at our 172 GeV sample is therefore smaller 
than that at 183 or 189 GeV, where the value for a s we obtain is larger. 
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determined using By/ at 183 and 189 GeV is somewhat smaller, as already observed at 
LEP-1 energies ||25|| , which may indicate that the higher order corrections to B\y are 
slightly different from those of the other event shapes. 

The correlations between the values of a s at the three different cm. energies are large. 
For example, the correlations between the theoretical uncertainties of the three values 
of a s quoted in equation || are all close to 100%. The correlation coefficients for the 
experimental uncertainties vary between 40% and 67%, and the ones for the hadronization 
corrections vary between 84% and 98%. 

In order to allow for the effects of these correlations, we construct a value of a s at the 
scale given by the luminosity weighted cm. energy of our three data samples, 186.7 GeV. 
We use O(a^) predictions to evolve all our values of a s to this scale. As a final result, 
using the weighted mean values of a s at all three energies, we obtain 

a s (187GeV) = 0.106 ± 0.001(stat.) ± 0.004(syst.). (4) 

When evolved to the scale of M z o, using 0(a^) calculations, the value of the weighted 
average for a s becomes a s (M z o) = 0.117 ± 0.005 (see also Table [7]). For comparison, 
our measurements at the M z o energy with a slightly different set of observables based on 



(9(a 2 )+NLLA QCD calculations yielded a s (M z0 ) = 0.120 ± 0.006 [|| . 

In order to compare average a s values using exactly the same theoretical predictions 
and observables as used in our previous publications, we compute the weighted average 
from our present fits to only 1 — T, M H , B T and B w . This is a s (187 GeV) = 0.104±0.005 
which corresponds to a s (M z o) = 0.115 ± 0.006 when evolved to the scale of M z o. Our 
previous analysis at ^/s = M z o gave a s (M z o) = 0.116 ± 0.006 from fits of the same 
predictions to the same observables. Our present determinations of a s are therefore 
consistent with our measurement at \/s = M z o. 

In Figure [8] we show our result ([$]) together with our other measurements of the strong 
coupling strength, as a function of the energy scale Q. The curve shows the 0(al) 
prediction for a s (Q) using a s (M z o) = 0.119 ± 0.004, the value of a s (M z o) as given in 
reference [[16|]. Note that this value of a s is obtained using many observables, from a 
wide variety of experimental environments. Our determinations of a s , obtained with the 



four event shapes mentioned above, all fall below the predictions of reference [46]. This 



indicates that the values obtained using these observables are somewhat shifted with 



respect to the value from reference [fig] . However, the consistency of our measurements 
between the LEP-1 and LEP-2 cm. energies is excellent. 

We obtained a relatively low value of a s using the 172 GeV sample alone for all event 
shape observables; approximately 1.8 standard deviations below expectations. 

In Figure ^ we show our results again, on a logarithmic scale, together with results 
obtained from the other experiments at lower cm. energies. 



5.3 Charged Multiplicity 

We measure the charged particle multiplicity distribution and derive several related quan- 
tities from it, in particular the mean charged multiplicity (n c h), the dispersion D = 
(,( n ch) ~ ( n ch) 2 )z, the ratio (n c b)/D, the normalised second moment C2 = (^ch)/( n ch) 2 
and the second binomial moment 0^2= (^ch(^ch — 1))/ (^ch) 2 - 
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To correct the observed charged particle multiplicity distribution, we first subtract 
the expected background from 4-fermion events as described in Section |4j. The result- 
ing distribution is corrected for experimental effects such as acceptance, resolution and 
secondary interactions in the detector with an unfolding matrix, as previously done in 
references [[j],[|7|]. This matrix is determined from the PYTHIA sample of fully simulated 
(Z°/7)* events. Biases introduced by the event selection, by radiative events passing the 
selection and by the fraction of particles with lifetimes shorter than 3 x 10~ 10 s that did 
not decay in the detector are corrected using a bin-by-bin multiplication method. The 
bin-by-bin corrections are typically smaller than 10-15%. 

The charged particle multiplicity distribution, corrected for experimental effects using 
an unfolding matrix, is tabulated in Table ||. The data are shown for our 189 GeV 
sample in Figure p!0|a. The COJETS model, as indicated in the figure, predicts too 
many high multiplicity events and clearly disagrees with the data. The PYTHIA and 
ARIADNE models predict mutually indistinguishable distributions which describe the 
data reasonably well. The HERWIG model gives a somewhat better description for low 
and intermediate values. 

We determine the mean values to be: 

(n ch )(172 GeV) = 25.77 ± 0.58(stat.) ± 0.88(syst.), 

(rich) (183 GeV) = 26.85 ± 0.27(stat.) ± 0.52(syst.), (5) 

(n ch )(189 GeV) = 26.95 ± 0.16(stat.) ± 0.51(syst.), 

The systematic errors are detailed in Table 9. The values of the dispersion D, the ratio 
(^ch) /D, the normalised second moment Ci and the second binomial moment 'R-i can be 
found in Table [5[ 

As a consistency check, the mean charged multiplicity is also computed by integrating 
the corrected rapidity distribution, and the fragmentation function, as determined in the 
next section. Both the rapidity distribution and the fragmentation function give results 
in good agreement with the one from the direct measurement presented in this section. 

Figure |10|b shows our measurements of (rich) together with results from OPAL at 



lower cm. energies. The data are compared to those of the other LEP experiments and 
to analytic QCD or Monte Carlo predictions. The predicted value of (n c h)(189 GeV) from 
PYTHIA is 27.6. This value changes by up to 0.4 when the hadronization parameters b 
and gq and the parton virtuality cutoff parameter Q are varied by ± 1 standard deviation 
about their tuned values ||15|| . The values from ARIADNE and HERWIG are 27.5 and 



27.2 respectively. The measured value is therefore about one standard deviations low 
compared to the PYTHIA and ARIADNE models. This trend is seen at all cm. energies 
of the LEP-2 data. However, this might be only a consequence of the high correlation 
between the systematic uncertainties of the OPAL measurements. COJETS predicts 
the charged multiplicity to be above 30, well above the measured value. The models 
have been tuned to agree with the OPAL data taken at the Z° peak with (n c h)(Mzo) = 
21.05 ± 0.01(stat.) ±0.20(syst.) @. 

The dispersion D of the data, D(189 GeV) = 8.45 ± 0.12(stat.) ± 0.34(syst.) is better 
described by the PYTHIA and ARIADNE Monte Carlo samples, which predict 8.65 and 
8.58 respectively. In contrast, the HERWIG model predicts _D(189 GeV) = 9.25, more 
than two standard deviations too large. 
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The dash-dotted curve in Figure |T0| b shows the NLLA QCD prediction for the 
energy evolution of the charged particle multiplicity, with parameters fitted to all available 
data points between 12 and 161 GeV [|1],|3],|^,|,0,|48],[5(|. The NLLA calculation predicts a 
mean charged multiplicity at 189 GeV of 27.6. Variations in the quark flavour composition 
at the different cm. energies do not significantly influence the results of the fit. 



5.4 Charged particle momentum spectra 

We measure the charged particle fragmentation function, 1/er • dcr c h/dx p , and the £ p dis- 
tribution, I/a • dcx c h/d£p, where x p = 2p/y/s, £ p = m(l/x p ) and p is the measured track 
momentum, using the same methods as in our previous publications |j],|2|]. We determine 
the rapidity distribution, y = | ln( ^T^jj )|, where p\\ is the momentum component parallel 
to the thrust axis and E is the energy of the particle. We also study the distribution of 
the 3-momentum components parallel, p±, and perpendicular, p° ut , to the event plane. 
This plane is defined by the eigenvectors of the momentum tensor associated with the two 
largest eigenvalues, as in reference 

The p™, p^ ut and y distributions for events with a cm. energy of yfs = 189 GeV are 



shown in Figure [11] and are tabulated in Tables [10], |TT] and We observe good agreement 
between PYTHIA, HERWIG and ARIADNE and the data for the p^ distributions. How- 
ever, the slopes of the distributions of the COJETS prediction are somewhat steeper 
than the data. In case of the p° ut distribution, not only COJETS but to lesser extent 
also PYTHIA, HERWIG and ARIADNE predict a slightly softer pf 1 spectrum than the 
spectrum observed in the data. In the y distribution, there is reasonable agreement of the 
PYTHIA, HERWIG and ARIADNE Monte Carlo models with the data. The COJETS 
Monte Carlo model significantly overestimates the production of charged particles with 
low values of y. 

The fragmentation function and the £ p distribution are shown in Figures [II] (d) 
and (a) together with the Monte Carlo predictions. Numerical values of these data 



are given in Tables [13] and [L4|. The spectrum of charged particles with large momentum 
fraction x p is well described by all Monte Carlo models. The shape of the £ p distribution 
is well described by PYTHIA, although their normalisation is somewhat higher, reflecting 
the difference between the predicted and the observed charged particle multiplicity (see 
Section |5.3| ). The COJETS Monte Carlo predicts too many particles in the region of the 
peak and at large values of £ p , where low momentum particles contribute. In the LLA 
approach, this is the region where soft gluon production is reduced as a consequence of 
destructive interference. 

Based on the concept of local parton hadron duality (LPHD) ]51] within the leading- 
log approximation of QCD calculations (LLA) |52|] one expects a Gaussian shape of the 
£ p distribution. The peak of the distribution is predicted to have an almost logarithmic 
variation with energy. The next to leading-log approximation of QCD calculations 
predicts the shape to be a Gaussian with higher moments (skewed Gaussian), and the 
same energy dependence of the peak as the LLA. We fit a skewed Gaussian of the form 
suggested in |53] to the region 2.0 < £ p < 6.2 of the three £ p distributions and determine 
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the positions of the peaks, £ 0) to be 



£ (172 GeV) = 4.031 ± 0.033(stat.) ± 0.041(syst.), 

£o(183 GeV) = 4.087 ± 0.014(stat.) ± 0.030(syst.), (6) 

^o( 189 GeV) = 4.124 ± 0.010(stat.) ±0.037(syst.). 

The systematic uncertainty takes into account experimental effects as described in Sec- 
tion [O] and the uncertainty due to the choice of the fit range. The fit range was reduced 
to 3.2 < £ p < 4.8 and increased to 1.6 < £ p < 6.6. The larger deviation from the result 
with the standard fit range was taken as systematic uncertainty and added in quadrature 
to the other effects. All systematic effects are summarised in Table 15. 

In Figure [E| (b) we show our measurements of £o together with measurements taken 
at various cm. energies 0^|,|||,|54|,|||, and the PYTHIA, HERWIG and CO JETS pre- 



dictions. The dash-dotted curve shows the prediction for the energy evolution by the 
modified leading-log approximation (MLLA) |56| with parameters fitted to the different 
energies, excluding the results of this paper. The values of the MLLA prediction extrap- 
olated to v^= 172, 183 and 189 GeV of 4.01 ±0.01, 4.05 ±0.01 and 4.07 ±0.01 are lower 
than our measurements. As noted in our previous publications [0,0], the fit is dominated 
by contributions from the data points with small errors at 29 and 35 GeV. The PYTHIA 
and HERWIG predictions are very similar to the MLLA fit, while the CO JETS prediction 
does not agree with the data. 

In Figure ^ (a) also the shape predicted by MLLA [J57]] is shown. The MLLA predic- 
tion is determined by two parameters, an effective QCD scale A e ff and a normalisation 
factor K{^fs) which within the framework of LPHD is expected to have a weak cm. en- 



ergy dependence. We fix A e jf to the value determined in f54}l , A e // = 0.253, and fit the 
formula to the measured charged particle momentum spectra in the region 2.8 < £ p < 4.8. 
The result of this fit to the ^/s =189 GeV data is shown in Figure [12] (a). The MLLA 
description and data show good agreement in the peak region but disagree at large £ p 
where kinematic effects become important and the perturbative QCD calculations are 



not valid 581. The values obtained for the normalisation factor K are: 



K(y/s = 172 GeV) = 1.143 ± 0.030(stat.) ± 0.056(syst.), 

K{\fs = 183 GeV) = 1.183 ± 0.010(stat.) ± 0.023(syst.), (7) 

K(yG = 189 GeV) = 1.164 ± 0.008(stat.) ± 0.030(syst.). 

To estimate the uncertainty due to the choice of the fit range, the fit range was reduced 
to 3.6 < £ p < 4.2 and increased to 2.4 < £ p < 5.2. the systematic uncertainty is shown in 



Table 16. In Figure |13|, these numbers are compared with results for the normalisation 
factors obtained in [54 for £ p distributions at lower cm. energies. The values from the 
new measurements are significantly below the value at yfs = M z o, K(M z o) = 1.28 ±0.01, 
thus following the observed trend of decreasing values for K with increasing cm. energy. 



A decrease of K values might be caused by higher order effects as discussed in . 

MLLA makes several predictions for relations between the different quantities that 
can be used to describe the centre of the £ p distribution, i.e. the position of the maximum 
£o, the mean value and the median value £ m which is defined as the value where 
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/Q m (da c h/d£p)d£ p is equal to /^(dcr c h/d£ p )d£ p . The difference £ — (£ p ) is predicted to be 
approximately 0.351 (0.355) for 3(5) active flavours, independent from the cm. energy 
|57| . The first non-leading corrections to this value cannot be predicted at the moment. 
However, this correction has been taken into account for the calculation of the ratio of 
differences, (£o — (£p))/(£m~ (£p)) which is predicted to be 3.5 (3.4) for 3(5) active flavours 
for the cm energies considered in this paper and assuming A e ff = 0.253. In our analysis, 
we measure 

£ - (£ P }(172 GeV) = 0.054 ± 0.023(stat.) ± 0.063(syst.), 

£o - (Cp)(183 GeV) = 0.012 ± 0.010(stat.) ± 0.064(syst.), (8) 

£o - (f P )(189 GeV) = -0.035 ± 0.007(stat.) ± 0.075(syst.) 

and 



(9) 



The systematic errors are detailed in Tables 17 and 18. The results strongly depend 
on the high £ p region of the £ p spectrum, i.e. on the low momentum region where the 
selection efficiency for tracks is low and the corrections large. In MLLA mass effects 
are neglected and the £ p distribution vanishes for values above £ p = log(E/A e ff). In 
56|.|6(J, a simple procedure was suggested on how to take these effects into account in 
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the MLLA predictions. Indeed, when we fitted these 'modified' MLLA predictions to the 
£ p distribution at yfs = 189 GeV, it turned out that the fitted values for £ , £ m and (£ p ) 
became numerically very close to each other such that the value for the difference £o — (£p) 
almost vanished, consistent with our measurement. This also indicates that the value and 
even the sign for the ratio (£o — (£ p ))/(£ m — (£ P )) is largely undetermined, and depends 
strongly on the assumptions made in the correction procedure to the MLLA distribution. 
A comparison with the measured values is therefore not conclusive. 



6 Summary and conclusions 

In this paper we have presented measurements of the properties of hadronic events pro- 
duced at LEP at centre-of-mass energies between 172 and 189 GeV. The 189 GeV data-set 
recorded by the OPAL detector constitutes approximately 70% of the total luminosity of 
the LEP-2 programme available by the end of 1998, and provides the most precise results 
at energies above the Z° resonance. 

We have determined the corrected distributions for event shape observables, for the 
charged particle multiplicity and for charged particle momentum spectra at centre-of-mass 
energies of ^7s =172, 183 and 189 GeV. The predictions of the PYTHIA, HERWIG and 
ARIADNE Monte Carlo models are found to be in general agreement with the measured 
distributions. While the COJETS Monte Carlo model describes most event shape distri- 
butions well, it overestimates the number of soft charged particles produced as observed 
in the rapidity and £ p = ln(l/x p ) distributions. 
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From a fit of 0(a^)+NLLA QCD predictions to five event shapes and the differential 
two jet rate, defined using the Durham jet finder, we have determined the strong coupling 
parameter a s (187 GeV) = 0.106 ± 0.001 (stat.) ± 0.004(syst.) at a luminosity weighted 
centre-of-mass energy of 186.7 GeV. When this is evolved to the Z° peak it is in excellent 
agreement with our previous measurement at s/s = M z o. 

The mean charged particle multiplicity has been determined to be (n c h)(189 GeV) = 
26.95 ± 0.16(stat.) ± 0.51(syst.), about one standard deviations below the NLLA QCD 
predictions and the PYTHIA and ARIADNE Monte Carlo predictions. 

The £ p distribution is not in perfect agreement with QCD MLLA calculations. The 
position of the peak was measured to be £ (189 GeV) = 4.124±0.010(stat.) ±0.037(syst.). 
This is approximately one and a half standard deviations larger than the expectation for 
the energy evolution given by a QCD MLLA extrapolation from the results of lower energy 
experiments. 

Our studies show that most of the features of hadronic events produced in e + e~ colli- 
sions at energies above the Z° mass are well described by QCD in the form of analytic or 
Monte Carlo predictions, or both. 
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172 GeV 


183 GeV 


189 GeV 


(cm. energy) 


(GeV) 


172.1 


182.7 


188.7 


Luminosity 


(Pb- 1 ) 


10.4 


57.2 


186.3 


Preselection 


1318 


6142 


18954 


'ISR-flt' Selection 




273 


1447 


4457 


4-fermion background 


% 


20.3 


26.6 


29.4 


Efficiency non-rad events 


% 


82.6 


83.3 


82.9 


Final selection 




228 


1098 


3277 


4-fermion background 


% 


7.3 


9.0 


10.1 


Efficiency non-rad events 


% 


79.4 


79.4 


78.5 



Table 1: Data samples at 172, 183 and 189 GeV. The luminosity and luminosity weighted 
mean cm. energies are given in the first two rows. The rows labelled as 'ISR-fit selection' 
and 'Final selection' correspond to the number of events passing these criteria. 
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T 


R(T) (172 GeV) 


R(T) (183 GeV) 


R(T) C189 GeV) 


n ^7nn n von 


0.62 ± 0.27 


± 0.31 


0.29 ± 0.11 


± 0.17 


n qi^ j_ n nn i n 1 a 

u.od ± u.uy ± u.i4 


U. /80-0.85U 


0.60 ± 0.28 


± 0.22 


0.75 ± 0.15 


± 0.16 


n 7o j_ n nn i n i a 

u. (o ± u.uy ± u.i4 


n Q^n n con 
0.850-0.880 


0.73 ± 0.45 


± 0.54 


1.45 ± 0.28 


± 0.44 


l.DO ± U.l I ± u.ZZ 


n con n m n 

0.880-0.910 


1.35 ± 0.52 


± 1.77 


1.62 ± 0.27 


± 0.68 


o oo _i_ nio _i_ ni^ 

2.22 ± U.lo ± U.lb 


n ni n n non 
0. 910-0. 93U 


3.01 ± 0.91 


1 O (\ 1 

± 2.91 


o on i r\ a a 

3.30 ± 0.44 


± 0.60 


o nn _i_ noc; j_ n ak 
o.Uy ± U.zo ± U.45 


u.you-u.you 


6.60 ± 1.25 


± 1.77 


5.27 ± 0.54 


± 0.64 


c i/i _i_nQi _i_ n A& 


950-0 960 

KJirJKJKJ KJ • ij KJ KJ 


4.47 ± 1.43 


± 3.20 


8.20 ± 0.95 


± 0.95 


7 02 =b 50 ± 88 


960-0 970 

v 7 . \7 \j v / \ i . t j t y i 


7.85 ± 1.90 


± 3.19 


10.0 ± 1.0 


± 1.4 


9 94 ± 59 ± 56 

/ • ■./ -L — 1 — V 7 • ■ / — 1 — V / . '.7 \J 


0.970-0.980 


18.5 ± 2.9 


± 6.8 


15.1 ± 1.3 


± 0.8 


15.6 ± 0.7 ± 1.2 


980-0 990 


22.3 ± 3.4 


± 3.7 


22.9 ± 1.7 


± 1.6 


22.0 ± 0.9 ± 1.6 


990-1 000 


11.6 ± 2.6 


± 3.8 


9.08 ± 1.07 


± 1.46 


8 73 ± 60 ± 63 


T 

- 1 major 


i?(T major ) (172 GeV) 


i?(T maj0 r) (183 GeV) 


major/ 1 v ioy vjjc v j 


n nnn n n/in 
U.UUU-U.U4(J 


1.68 ± 0.57 


± 0.62 


0.92 ± 0.19 


± 0.16 


n qq i n 11 i n 1 a 
U.oo ± U.ll ± U.14 


n n/i n n ncn 
U.U4U-U.U8U 


7.33 ± 0.99 


± 0.68 


7.17 ± 0.47 


± 0.34 


{? CZR _l n O/^ _l n QK 
D.DD ± U.ZD ± U.OO 


n ncn n i on 
U.UoU-U.lzU 


4.87 ± 0.72 


± 0.18 


4.96 ± 0.36 


± 0.27 


k nn _l n 01 _l n oo 
O.UU ± U.Zl ± [J.22 


o 1 90 n i fin 

U. ±ZrU _ U. 1UU 


3.23 ± 0.60 


± 1.09 


3.22 ± 0.29 


± 0.60 


Q1Q_Lni7_LnQQ 

o.iy zc u.i/ zc u.oo 


160-0 220 


2.03 ± 0.40 


± 0.28 


2.32 ± 0.20 


± 0.19 


2 36 ± 12 ± 20 


0.220-0.300 


1.07 ± 0.27 


± 0.38 


1.18 ± 0.14 


± 0.24 


1.26 ± 0.08 ± 0.14 


300-0 400 

KJ > ' f \ 7 \ / V / ■ A. y 7 \ 7 


0.49 ± 0.21 


± 0.21 


0.84 ± 0.12 


± 0.21 


73 ± 07 ± 05 

U . 1 '.7 — 1 — \_7 . KJ 1 — 1 — KJ . KJ .7 


400-0 500 

v7 • JI V7 v7 V7«t7V7W 


0.29 ± 0.18 


± 0.13 


0.23 ± 0.08 


± 0.12 


35 ± 06 ± 12 

v / . 1 7 '.J — 1 — KJ • KJ KJ — 1 — V7 • -L — J 


500-0 600 

v / . ' 7 v / v y v / > \7 y 7 \ 7 


0.33 ± 0.17 


± 0.24 


0.033 ± 0.075 ± 0.124 


18 ± 06 ± 11 


T ■ 

± minor 


R(T minor ) (172 GeV) 


i?(T minor ) (183 GeV) 


R(T ■ ) C189 GeVl 


n nnn n non 
0.000-0.020 


0.59 ± 0.43 


± 0.54 


0.087 ± 0.067 ± 0.129 


n oo _L n n^ _i_ n 1 o 
U.22 ± U.Ub ± U.12 


n non n n /i n 
0.02U-0.04U 


10.6 ± 1.8 


± 1.5 


9.49 ± 0.80 


± 0.82 


q no _i_ n /ii _i_ i 9c 
8.03 ± 0.41 ± 1.38 


n n/i n n n/^n 
U.U4U-U.U0U 


16.8 ± 2.1 


± 5.6 


15.7 ± 1.0 


± 1.5 


14. / ± U.O ± 1. ( 


n n^n n non 
0.060-0.080 


7.73 ± 1.31 


± 2.77 


10.6 ± 0.8 


± 1.8 


1U.1 ± U.4 ± 1.6 


n ncn n 1 nn 
U.UoU-U.lUU 


6.24 ± 1.16 


± 2.16 


5.86 ± 0.55 


± 1.26 


k ka -L- n Q1 _i_ n F^n 
0.04 ± U.ol ± U.OU 


n mn n 19(1 


4.07 ± 0.99 


± 2.07 


3.46 ± 0.44 


± 0.72 


9 nc 7 n OA -4- D 94 
l.\7o in u.z^ m u.z^i 


120-0 140 


1.52 ± 0.66 


± 1.45 


2.69 ± 0.44 


± 0.86 


2 06 ± 22 ± 32 


140-0 160 

V 7 . -L A V 7 V / > J. V / V / 


1.16 ± 0.67 


± 0.97 


1.24 ± 0.33 


± 0.26 


1.55 ± 0.21 ± 0.36 


160-0 200 

v / . _l y.j y 7 v / ■ _ v / v / 


0.48 ± 0.39 


± 0.77 


0.69 ± 0.23 


± 0.66 


0.69 ± 0.14 ± 0.44 


200-0 240 


0.16 ± 0.32 


± 0.76 


0.20 ± 0.22 


± 0.55 


12 ± 15 ± 75 

KJ • J_ j-J — 1 — V 7 . 1-7 — 1 — KJ • 1 KJ 


240-0 300 


0.30 ± 0.48 


± 0.86 


-.475 ± 0.617 ± 1.540 


15 ± 18 ± 46 

V7 . X t7 — 1 — V7 • KJ — 1 — V ^ • i \.7 


A 


R(A) (172 GeV) 


R(A) (183 GeV) 


R(A) f IRQ GeV) 


n nnn n nnc; 
U.UUU-U.UUo 


119. ± 10. 


± 34. 


131. ± 5. 


± 5. 


1 on i o i n 

Izy. ± o. ± y. 


n oo^ o m n 


33.2 ± 5.9 


± 9.8 


34.2 ± 3.1 


± 3.3 


on q _l i 7 j_ i o 
ou.y m i.i ziz x.o 


0.010-0.015 


14.2 ± 4.0 


± 9.4 


12.6 ± 2.0 


± 4.4 


12.7 ± 1.2 ± 0.8 


0.015-0.025 


5.88 ± 2.15 


± 2.02 


5.56 ± 1.07 


± 2.28 


4.87 ± 0.62 ± 1.40 


0.025-0.040 


0.80 ± 0.91 


± 1.42 


2.55 ± 0.72 


± 1.26 


2.17 ± 0.44 ± 0.58 


0.040-0.070 


-.442 ± 0.336 ± 0.400 


0.58 ± 0.44 


± 1.42 


0.59 ± 0.25 ± 0.43 


0.070-0.100 


0.14 ± 0.44 


± 7.30 


-.419 ± 0.620 ± 1.281 


0.14 ± 0.32 ± 0.66 



Table 2: Results R(y) — 1/a • da/dy at y/s = 172, 183 and 189 GeV for the event shape 
observables y: thrust T, thrust major T ma j or , thrust minor T minor and aplanarity A. The 
first error is statistical, the second systematic. 
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o 



R(0) (172 GeV) 



R(0) (183 GeV) 



R(0) (189 GeV) 



0.000- 
0.050- 
0.100- 
0.150- 
0.200- 
0.250- 
0.300- 
0.400- 



0.050 
0.100 
0.150 
0.200 
0.250 
0.300 
0.400 
0.500 



10.6 ± 
4.63 ± 
1.93 ± 
0.91 ± 
0.46 ± 
0.51 ± 
0.28 ± 
0.21 ± 



1.0 

0.68 

0.47 

0.34 

0.30 

0.29 

0.15 

0.11 



± 0.5 
± 0.71 
± 0.23 
± 0.23 
± 0.37 
± 0.19 
± 0.21 
± 0.15 



10.2 ± 
4.30 ± 
2.18 ± 
1.32 ± 
0.78 ± 
0.72 ± 
0.21 ± 
0.031 ± 



0.5 

0.31 

0.23 

0.20 

0.16 

0.16 

0.07 



± 0.4 
± 0.51 
± 0.25 
± 0.57 
± 0.47 
± 0.19 
± 0.09 



0.034 ± 0.062 



9.97 ± 
4.15 ± 

2.33 ± 

1.34 ± 
0.82 ± 
0.60 ± 
0.29 ± 

0.087 ± 



0.27 
0.18 
0.14 
0.12 
0.10 
0.09 
0.05 



± 0.44 
± 0.38 
± 0.11 
± 0.23 
± 0.21 
± 0.18 
± 0.10 



0.028 ± 0.057 



C 



R(C) (172 GeV) 



R(C) (183 GeV) 



R{C) (189 GeV) 



0.000-0.050 
0.050-0.080 
0.080-0.110 
0.110-0.140 
0.140-0.180 
0.180-0.220 
0.220-0.300 
0.300-0.400 
0.400-0.500 
0.500-0.600 
0.600-0.750 
0.750-1.000 



2.76 ± 

4.65 ± 
4.68 ± 

3.66 ± 
2.12 ± 
1.31 ± 
1.62 ± 
0.82 ± 
0.28 ± 
0.28 ± 
0.45 ± 

-.002 ± 



0.59 
0.92 
0.86 
0.73 
0.48 
0.39 
0.31 
0.21 
0.14 
0.15 
0.16 



± 1.12 
± 2.48 
± 2.48 
± 0.86 
± 0.86 
± 0.48 
± 0.35 
± 0.21 
± 0.38 
± 0.26 
± 0.24 



0.057 ± 0.153 



1.78 ± 
5.70 ± 
3.85 ± 
3.31 ± 
2.28 ± 
1.77 ± 
1.46 ± 
0.80 ± 
0.60 ± 
0.44 ± 
0.26 ± 
0.093 ± 



0.22 
0.49 
0.37 
0.34 
0.24 
0.21 
0.14 
0.09 
0.09 
0.08 
0.07 



± 0.33 
± 0.38 
± 0.60 
± 0.43 
± 0.33 
± 0.25 
± 0.13 
± 0.10 
± 0.20 
± 0.19 
± 0.09 



0.058 ± 0.102 



1.82 ± 
4.97 ± 
4.21 ± 
3.14 ± 
2.45 ± 
1.80 ± 
1.27 ± 
0.86 ± 
0.66 ± 
0.42 ± 
0.35 ± 
0.082 ± 



0.13 
0.26 
0.22 
0.19 
0.15 
0.13 
0.08 
0.06 
0.05 
0.05 
0.05 



± 0.17 
± 0.19 
± 0.67 
± 0.47 
± 0.24 
± 0.28 
± 0.13 
± 0.08 
± 0.08 
± 0.14 
± 0.07 



0.035 ± 0.060 



R(M H ) (172 GeV) 



R(M H ) (183 GeV) 



R(M H ) (189 GeV) 



0.060-0.075 
0.075-0.090 
0.090-0.110 
0.110-0.140 
0.140-0.170 
0.170-0.200 
0.200-0.250 
0.250-0.300 
0.300-0.350 
0.350-0.450 
0.450-0.600 



1.57 ± 
5.52 ± 
7.15 ± 
5.84 ± 
4.27 ± 
4.03 ± 
2.70 ± 
0.88 ± 
0.72 ± 
0.56 ± 
0.34 ± 



0.53 
1.26 
1.37 
0.99 
0.84 
0.84 
0.54 
0.36 
0.34 
0.22 
0.13 



± 1.32 
± 1.42 
± 2.60 
± 1.69 
± 1.44 
± 1.18 
± 0.50 
± 0.64 
± 0.21 
± 0.09 
± 0.19 



1.20 ± 
4.95 ± 
6.31 ± 
5.62 ± 
5.06 ± 
3.88 ± 
2.71 ± 
1.56 ± 
1.08 ± 
0.57 ± 
0.11 ± 



0.23 
0.57 
0.61 
0.48 
0.44 
0.39 
0.25 
0.21 
0.19 
0.11 
0.05 



± 0.47 
± 0.98 
± 0.88 
± 0.63 
± 0.80 
± 0.66 
± 0.42 
± 0.14 
± 0.32 
± 0.04 
± 0.07 



1.67 ± 
4.35 ± 
6.07 ± 
6.42 ± 
4.77 ± 
3.22 ± 
2.45 ± 
1.66 ± 
1.37 ± 
0.57 ± 
0.14 ± 



0.17 
0.31 
0.34 
0.29 
0.25 
0.20 
0.14 
0.13 
0.13 
0.07 
0.03 



± 0.23 
± 0.48 
± 0.56 
± 0.39 
± 0.35 
± 0.32 
± 0.15 
± 0.11 
± 0.19 
± 0.14 
± 0.10 



R(S) (172 GeV) 



R(S) (183 GeV) 



R(S) (189 GeV) 



0.000-0.020 
0.020-0.040 
0.040-0.060 
0.060-0.120 
0.120-0.200 
0.200-0.300 
0.300-0.500 
0.500-0.700 



26.7 ± 
8.19 ± 
3.63 ± 
1.43 ± 
0.60 ± 
0.24 ± 
0.17 ± 
0.19 ± 



2.5 

1.43 

0.95 

0.39 

0.25 

0.19 

0.11 

0.12 



± 2.0 
± 1.60 
± 0.60 
± 0.26 
± 0.41 
± 0.22 
± 0.16 
± 0.25 



24.7 ± 
7.70 ± 
4.18 ± 
2.04 ± 
0.65 ± 
0.44 ± 
0.093 ± 
0.16 ± 



± 2.0 
± 1.16 
± 0.71 
± 0.16 
± 0.17 
± 0.22 
0.060 ± 0.093 
0.08 ± 0.13 



1.1 

0.65 
0.51 
0.21 
0.12 
0.11 



25.0 ± 
7.46 ± 
4.15 ± 
1.90 ± 
0.79 ± 
0.37 ± 
0.22 ± 
0.055 ± 



0.7 

0.38 

0.29 

0.12 

0.08 

0.07 

0.05 



± 0.8 
± 0.69 
± 0.49 
± 0.12 
± 0.18 
± 0.04 
± 0.07 



0.030 ± 0.046 



Table 3: Results R(y) = l/a- da/dy &t y/s = 172, 183 and 189 GeV for the event shape 
observables y: oblateness O, C-parameter C, heavy jet mass M H and sphericity S. The 
first error is statistical, the second systematic. 
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R(B T ) (172 GeV) 


R(B T 


) (183 GeV) 


R(B T 


) (189 GeV) 


0.000-0.030 


3.96 


± 


0.98 


± 1.68 


1.40 ± 


0.26 


± 0.60 


1.89 


± 


0.17 


± 


0.36 


0.030-0.040 


11.3 


± 


2.6 


± 2.6 


13.4 ± 


1.3 


± 2.0 


9.84 


± 


0.64 


± 


0.56 


0.040-0.050 


8.75 


± 


2.09 


± 7.73 


12.9 ± 


1.3 


± 2.7 


11.9 


± 


0.7 


± 


1.7 


0.050-0.060 


11.7 


± 


2.3 


± 9.7 


8.57 ± 


0.91 


± 1.54 


9.59 


± 


0.57 


± 


0.46 


0.060-0.075 


8.82 


± 


1.59 


± 2.61 


8.69 ± 


0.78 


± 0.96 


8.00 


± 


0.42 


± 


0.53 


0.075-0.090 


4.46 


± 


1.14 


± 3.65 


5.20 ± 


0.58 


± 0.99 


6.33 


± 


0.38 


± 


0.55 


0.090-0.110 


5.53 


± 


1.12 


± 1.06 


5.67 ± 


0.56 


± 0.91 


4.37 


± 


0.27 


± 


0.61 


0.110-0.130 


3.19 


± 


0.85 


± 1.47 


3.11 ± 


0.40 


± 0.64 


3.88 


± 


0.26 


± 


0.39 


0.130-0.160 


2.09 


± 


0.59 


± 0.92 


2.64 ± 


0.32 


± 0.29 


2.55 


± 


0.18 


± 


0.24 


0.160-0.200 


1.21 


± 


0.46 


± 0.66 


1.64 ± 


0.24 


± 0.29 


1.75 


± 


0.15 


± 


0.09 


0.200-0.250 


0.49 


± 


0.35 


± 0.45 


0.93 ± 


0.21 


± 0.36 


0.76 


± 


0.13 


± 


0.47 


0.250-0.300 


1.11 


± 


0.46 


± 0.69 


0.27 ± 


0.20 


± 0.38 


0.66 


± 


0.16 


± 


0.25 


0.300-0.350 


-.002 


± 


0.345 ± 0.866 


0.40 ± 


0.30 


± 0.52 


0.058 


± 


0.149 ± 


0.219 


B w 


R{B W 


) (172 GeV) 


R{B W 


) (183 GeV) 


R{B W 


) (189 GeV) 


0.000-0.020 


8.91 


± 


1.75 


± 3.96 


4.81 ± 


0.59 


± 0.93 


4.61 


± 


0.34 


± 


0.91 


0.020-0.030 


14.5 


± 


2.7 


± 4.0 


18.1 ± 


1.5 


± 1.0 


16.8 


± 


0.8 


± 


1.8 


0.030-0.040 


13.7 


± 


2.4 


± 4.5 


13.1 ± 


1.2 


± 1.0 


13.7 


± 


0.7 


± 


0.5 


0.040-0.050 


8.48 


± 


1.89 


± 3.88 


10.5 ± 


1.0 


± 0.4 


10.3 


± 


0.6 


± 


0.8 


0.050-0.065 


7.36 


± 


1.49 


± 2.67 


8.23 ± 


0.76 


± 0.86 


7.89 


± 


0.43 


± 


0.42 


0.065-0.080 


6.41 


± 


1.41 


± 2.33 


5.16 ± 


0.61 


± 1.06 


5.44 


± 


0.37 


± 


0.35 


0.080-0.100 


3.88 


± 


1.01 


± 1.43 


4.51 ± 


0.51 


± 0.51 


4.30 


± 


0.30 


± 


0.32 


0.100-0.150 


1.57 


± 


0.46 


± 0.52 


2.30 ± 


0.26 


± 0.38 


2.66 


± 


0.17 


± 


0.27 


0.150-0.200 


1.02 


± 


0.43 


± 0.31 


1.25 ± 


0.24 


± 0.28 


1.03 


± 


0.13 


± 


0.24 


0.200-0.250 


0.77 


± 


0.35 


± 0.39 


0.33 ± 


0.14 


± 0.21 


0.53 


± 


0.11 


± 


0.23 


0.250-0.300 


0.048 


± 


0.099 ± 0.331 


0.014 ± 


0.058 ± 0.195 


0.048 


± 


0.038 ± 


0.074 


v& 


R{y 


^ (172 GeV) 


R(v& 


) (183 GeV) 


R(y 


23 


) (189 GeV) 


0.00030-0.00075 


399. 


± 


70. 


± 139. 


335. ± 


31. 


± 56. 


348. 


± 


18. 


± 


56. 


0.00075-0.00130 


202. 


± 


42. 


± 40. 


215. ± 


21. 


± 25. 


163. 


± 


10. 


± 


31. 


0.00130-0.00230 


84.5 


± 


19.7 


± 41.5 


101. ± 


10. 


± 9. 


111. 


± 


6. 


± 


10. 


0.00230-0.00400 


56.4 


± 


12.1 


± 24.6 


56.3 ± 


5.9 


± 9.8 


62.5 


± 


3.6 


± 


7.5 


0.00400-0.00700 


34.0 


± 


7.2 


± 5.6 


29.7 ± 


3.2 


± 5.1 


29.7 


± 


1.9 


± 


6.4 


0.00700-0.01200 


16.8 


± 


4.2 


± 6.0 


18.7 ± 


2.1 


± 5.4 


16.6 


± 


1.1 


± 


2.6 


0.01200-0.02300 


8.20 


± 


2.10 


± 1.99 


7.77 ± 


0.96 


± 1.50 


9.01 


± 


0.61 


± 


1.14 


0.02300-0.04000 


2.56 


± 


0.98 


± 1.74 


4.67 ± 


0.60 


± 1.06 


4.21 


± 


0.35 


± 


0.96 


0.04000-0.07000 


2.46 


± 


0.73 


± 1.47 


1.74 ± 


0.30 


± 0.62 


1.98 


± 


0.19 


± 


0.25 


0.07000-0.13000 


0.19 


± 


0.24 


± 1.02 


0.72 ± 


0.16 


± 0.28 


0.76 


± 


0.10 


± 


0.08 


0.13000-0.23500 


0.45 


± 


0.19 


± 0.28 


0.19 ± 


0.07 


± 0.08 


0.27 


± 


0.05 


± 


0.10 


0.23500-0.40000 


0.036 


± 


0.048 ± 0.074 


0.050 ± 


0.054 ± 0.111 


0.014 


± 


0.018 ± 


0.037 



Table 4: Results R(y) — 1/a ■ da/dy at y/s = 172, 183 and 189 GeV for the event shape 
observables y: total jet broadening B T , wide jet broadening B w and the transition value 
between 2- and 3-jets yf? 3 - The first error is statistical, the second systematic. 
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172 GeV 


1 — T 


M H 


Bt 


B\y 


C 


3/23 


Mean 


fit range 


0.05-0.30 


0.17-0.45 


0.075-0.25 


0.05-0.20 


0.18-0.60 


0.0023-0.130 


_ 


a s 


0.0999 


0.0912 


0.0908 


0.0897 


0.0910 


0.0976 


0.0919 


X 2 /clof 


6.7/5 


5.0/4 


3.3/5 


2.7/4 


5.9/4 


7.9/6 


- 


Stat 


±0.0094 


±0.0083 


±0.0076 


±0.0061 


±0.0095 


±0.0073 


±0.0064 


Exp 


±0.0127 


±0.0084 


±0.0085 


±0.0061 


±0.0100 


±0.0203 


±0.0074 


Hadr 


±0.0019 


±0.0016 


±0.0032 


±0.0014 


±0.0023 


±0.0028 


±0.0015 


x M = 0.5 


-0.0030 


-0.0015 


-0.0029 


-0.0018 


-0.0024 


±0.0000 


-0.0020 


•Eft 2 


±0.0038 


±0.0023 


±0.0036 


±0.0025 


±0.0031 


±0.0016 


±0.0028 


Syst 


±0.0133 


±0.0087 


±0.0096 


±0.0066 


±0.0106 


±0.0205 


±0.0079 


Total 


±0.0163 


±0.0121 


±0.0123 


±0.0090 


±0.0142 


±0.0217 


±0.0102 


183 GeV 


1 — T 


M H 


Bt 


B\y 


C 




Mean 


fit range 


0.05-0.30 


0.17-0.45 


0.075-0.25 


0.05-0.20 


0.18-0.60 


0.0023-0.130 


_ 


a s 


0.1070 


0.1074 


0.1050 


0.1001 


0.1088 


0.1073 


0.1056 


X 2 /doi 


3.6/5 


2.9/4 


6.8/5 


2.8/4 


1.8/4 


3.0/6 


- 


Stat 


±0.0032 


±0.0036 


±0.0029 


±0.0025 


±0.0039 


±0.0032 


±0.0027 


Exp 


±0.0033 


±0.0032 


±0.0056 


±0.0027 


±0.0042 


±0.0024 


±0.0019 


Hadr 


±0.0015 


±0.0017 


±0.0021 


±0.0008 


±0.0016 


±0.0016 


±0.0010 


x M = 0.5 


-0.0038 


-0.0028 


-0.0045 


-0.0030 


-0.0043 


-0.0004 


-0.0026 


x^ = 2 


±0.0048 


±0.0040 


±0.0055 


±0.0038 


±0.0054 


±0.0024 


±0.0039 


Syst 


±0.0056 


±0.0050 


±0.0078 


±0.0044 


±0.0066 


±0.0032 


±0.0039 


Total 


±0.0065 


±0.0062 


±0.0083 


±0.0051 


±0.0077 


±0.0045 


±0.0047 


189 GeV 


1 — T 


M H 


Bx 


B\y 


C 




Mean 


fit range 


0.05-0.30 


0.17-0.45 


0.075-0.25 


0.05-0.20 


0.18-0.60 


0.0023-0.130 


_ 


a s 


0.1122 


0.1053 


0.1064 


0.1011 


0.1086 


0.1095 


0.1067 


X 2 /dof 


2.0/5 


5.1/4 


7.6/5 


5.1/4 


1.7/4 


6.4/6 


- 


Stat 


±0.0016 


±0.0017 


±0.0014 


±0.0013 


±0.0018 


±0.0017 


±0.0013 


Exp 


±0.0021 


±0.0027 


±0.0041 


±0.0020 


±0.0019 


±0.0036 


±0.0021 


Hadr 


±0.0013 


±0.0018 


±0.0020 


±0.0008 


±0.0014 


±0.0017 


±0.0008 


x M = 0.5 


-0.0045 


-0.0029 


-0.0047 


-0.0031 


-0.0044 


-0.0004 


-0.0030 


x^ = 2 


±0.0056 


±0.0040 


±0.0056 


±0.0040 


±0.0054 


±0.0026 


±0.0043 


Syst 


±0.0056 


±0.0047 


±0.0069 


±0.0041 


±0.0054 


±0.0043 


±0.0043 


Total 


±0.0058 


±0.0050 


±0.0070 


±0.0043 


±0.0057 


±0.0046 


±0.0045 



Table 5: Values of a s (172 GeV), a s (183 GeV) and a s (189 GeV) derived using the 
0(a 2 s )+~NLLA QCD calculations with x^ — 1 and the ln(i?)-matching scheme, fit ranges 
and x 2 /d.o.f. values for each of the six event shape observables. In addition, the statistical, 
and various systematic uncertainties are given. The sign indicates the direction in which 
a s changes with respect to the standard analysis. 
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1 79 CeV 

11— vjrc v 


1 PrPV 

lOO Ore V 


1 HQ ripV 


1 87 CIpV 

XO 1 VJTC V 


at s {Q) 


n nn 1 n 

u.uyiy 


n incc 
U.lUoo 


n 1 fifcT 
U.lUo ( 


n 1 n/^n 
U.lUbU 


Statistical error 


±U.UU04 


i n nno7 

ztU.UUZ i 


in nni q 


in nni q 


Iracks+Glusters 


n 1 t 

—0.0017 


+0.0004 


i n nnn r 

+0.0005 


i n nnn /i 

+0.0004 


tracks Unly 


i n nn i n 

+(J.UU19 


i n nnno 

+U.UUU8 


i n nn 1 n 

+U.UU19 


i n nm 7 

+U.UU17 


coster < U.7 


n nn /i n 

— U.UU49 


i n nnno 

+U.UUU2 


n nnno 

— U.UUU3 


n nnno 

— (J.UUU3 


Alternative v s 


i a a a o K 

+0.0035 


i a a a -i c\ 

+0.0012 


+0.0007 


i A A A A A 

+0.0009 


Wqcd 


+0.0024 


i a a a -i c\ 

+0.0012 


■ a a a a o 

+0.0003 


i A AAA i~\ 

+0.0002 


Background ±5/o 


i a a a a ,1 

+0.0004 


+0.0002 


i A AAA i~\ 

+0.0002 


i A AAA i~\ 

+0.0002 


Experimental syst. 


i a a a a 

±0.0074 


i a a a a 

±0.0019 


1 A A A (~\ -1 

±0.0021 


1 A A A (~\ A 

±0.0020 


6-1 s.d. 


-0.0001 


-0.0001 


-0.0001 


-0.0001 


6 + 1 s.d. 


+0.0001 


+0.0001 


+0.0001 


+0.0001 


-i 1 

a q — 1 s.d. 


+0.0002 


i a a a a -i 

+0.0001 


i A A A A -1 

+0.0001 


i A A A A -1 

+0.0001 


CTq + 1 S.d. 


—0.0002 


—0.0001 


A AAA -1 

—0.0001 


A AAA -1 

—0.0001 




_i_n nnm 


_i_n nnnzL 


_i_n nnn^ 
-rU.uuuo 


in firing 
-rU.uuuo 


ndsc onlv 


+0 0010 


+0 0003 


+0 0002 


+0 0003 


HERWIG 5.9 


-0.0006 


-0.0003 


-0.0002 


-0.0003 


ARIADNE 4.08 


+0.0009 


+0.0007 


+0.0007 


+0.0007 


Total hadronisation 


±0.0015 


±0.0010 


±0.0008 


±0.0009 


Xf, = 0.5 


-0.0020 


-0.0026 


-0.0030 


-0.0029 




+0.0028 


+0.0039 


+0.0043 


+0.0041 


Tot. syst. error 


±0.0079 


±0.0039 


±0.0043 


±0.0041 



Table 6: Weighted mean values of a s , derived from fits to six event shapes using the 
(9(a^)+NLLA QCD calculations with x M = 1 and the In (R) -matching scheme, for cm. 
energies at 172, 183 and 189 GeV. Full detailed breakdown of the systematic uncertainties 
are listed. The sign indicates the direction in which a s changes with respect to the 
standard analysis. The last column correspond to the values of a s , constructed at the 
luminosity weighted mean cm. energy. 





172 GeV 


183 GeV 


189 GeV 


187 GeV 


a s (M z o) 


0.0992 


0.1165 


0.1183 


0.1173 


Statistical error 
Total error 


±0.0072 
±0.0120 


±0.0033 
±0.0057 


±0.0016 
±0.0056 


±0.0015 
±0.0053 



Table 7: Weighted mean values of a s , as determined in Table ||, evolved to M z o. 
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P(n ch ) 172 GeV 


[%] 


P(n ch ) 183 GeV 


[%] 


P(n ch ) 189 GeV 


[%] 


8. 


94 


4- 


1 Q 


4_ 

31 


.28 


nq 


4- 


.uo 


4- 

31 


.08 
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. 1 1 


4_ 
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31 


.05 


10. 
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.uo 


-L 
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4- 
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4_ 


.13 
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4_ 
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4- 
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o .o j. 


i 


83 
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4- 
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4- 


36 


4- 


.49 


5 1 3 
o . ±o 


4- 


91 

. _■ J. 


4- 
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1 01 


4- 
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1 .1)1 


4- 
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4- 
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1 . ± O 


4- 


96 


4- 
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20. 
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4- 


1 1 3 


4- 


1.15 
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y.o i 


4- 


^9 

.OA 


4- 


.72 


R fi7 


4- 


9S 


4- 


.41 


22. 




4- 


1 09 


4- 


.98 


1 40 


4- 


53 

. OO 


4- 


.56 


9 46 


4- 


98 

. z*o 


4- 


.43 


24. 


i n ^9 


-U 


i 


4- 


.78 


10 3^ 


4- 
i 


. OO 


4- 


.47 


Q 79 


4- 
i 


9Q 


4- 


.42 


26. 


q i n 

y . ±u 


4- 


Q9 
. y^ 


4- 


.86 


Q fi8 
y.uo 


4- 


4Q 


4- 


.51 


Q A A 


4- 


98 


4- 


.36 


28. 


Q 1 3 
y . io 


4- 


qq 

. y y 


4- 


.85 




4- 


A^ 


4- 


.60 


R 77 


4- 


97 


4- 


.26 


30. 


7 73 


4- 


88 

. OO 


4- 


.76 


7 30 

1 . ou 


4- 


43 
.'-to 


4- 


.66 


7 87 

1 . 1 


4- 


9^ 

. z*o 


4- 


.25 


32. 


39 


-U 

1 


7q 
. t y 


4- 


.75 


on 

u.uu 


4- 
i 


.o t 


4- 


.55 


70 
u. / u 


4- 


9^ 

.zO 


4- 


.25 


34. 


^ 01 


4- 


69 


4- 


.83 


L ±. -L 


4- 


34 

. O 1 


4- 


.37 


O.O 1. 


4- 


99 


4- 


.26 


36. 


o.uy 


4- 


^7 


4- 


.97 


o.ou 


4- 


39 

.Oz 


4- 


.27 


4 ^4 


4- 


1 Q 


4- 

1 


.26 


38. 


9 49 


4- 


44 

. L ± L ± 


-U 

1 


.86 


3 1 3 

O. J.O 


4- 


31 


4- 


.25 


o.ou 


4- 


1 7 


4- 


.30 


40. 


1.71 


± 


.37 


± 


.68 


2.61 


± 


.30 


± 


.26 


2.52 




.15 


± 


.31 


42. 


1.19 


± 


.31 


± 


.45 


2.11 


± 


.27 


± 


.21 


1.88 


± 


.13 


± 


.27 


44. 


.55 


± 


.16 


± 


.26 


1.61 


± 


.26 


± 


.19 


1.35 


± 


.12 


± 


.24 


46. 


.36 


± 


.16 


± 


.22 


1.15 


± 


.27 


± 


.14 


1.00 


± 


.12 


± 


.21 


48. 


.30 


± 


.25 


± 


.26 


.88 


± 


.21 


± 


.16 


.72 


± 


.11 


± 


.21 


50. 


.35 


± 


.38 


± 


.30 


.56 


± 


.22 


± 


.21 


.50 


± 


.10 


± 


.18 


52. 


.07 


± 


.07 


± 


.29 


.32 


± 


.22 


± 


.24 


.34 


± 


.09 


± 


.16 


54. 


.01 


± 


.02 


± 


.28 


.19 


± 


.11 


± 


.19 


.21 


± 


.09 


± 


.12 


56. 


.02 


± 


.04 


± 


.31 


.14 


± 


.12 


± 


.14 


.12 


± 


.06 


± 


.09 


58. 


.26 


± 


.39 


± 


.25 


.00 


± 


.00 


± 


.12 


.06 


± 


.04 


± 


.11 


60. 


.01 


± 


.03 


± 


.72 


.03 


± 


.04 


± 


.13 


.09 


± 


.10 


± 


.09 


62. 


.01 


± 


.02 


± 


.28 


.00 


± 


.00 


± 


.13 


.02 


± 


.06 


± 


.09 


(n c h) 


25.77 ± 


0.58 ± 


0.88 


26.85 ± 


0.27 ± 


0.52 


26.95 ± 


0.16 ± 


0.51 


D 


8.38 ± 


0.63 ± 


1.29 


8.36 ± 


0.20 ± 


0.46 


8.45 ± 


0.12 ± 


0.34 


M/D 


3.08 ± 


0.21 ± 


0.45 


3.21 ± 


0.07 ± 


0.14 


3.19 ± 


0.04 ± 


0.09 


c 2 


1.106 ± 0.014 ± 0.030 


1.097 ± 0.004 ± 0.017 


1.098 ± 0.003 ± 0.005 


K 2 


1.067 ± 0.014 ± 0.031 


1.060 ± 0.004 ± 0.017 


1.061 ± 0.003 ± 0.006 



Table 8: Charged particle multiplicity in percent, measured at yfs = 172, 183 and 
189 GeV. The first error is statistical, the second systematic. 



29 





{n c h) 


172 GeV 


183 GeV 


189 GeV 


standard result 


25.77 


26.85 


26.95 


statistical error 


0.58 


0.27 


0.16 


Alternative \fp 


-0.15 


+0.01 


-0.01 


Background ±5% 


+0.01 


+0.02 


-0.02 


\d \ < 5 cm 


+0.29 


+0.05 


+0.05 


\zo\ < 10 cm 


+0.33 


+0.35 


+0.34 


N hits >S0 


-0.07 


-0.08 


-0.08 


| cos6> T | < 0.7 


-0.05 


+0.12 


+0.20 


HERWIG 


<0.01 


-0.35 


-0.20 


Wqcd 


+0.74 


-0.05 


+0.23 


Tot. syst. error 


0.88 


0.52 


0.51 



Table 9: Results with statistical and systematic uncertainties for the mean charged 
multiplicity (n ch ) at ^ = 172, 183 and 189 GeV. 



R( P f) (172 GeV) 



R&l) (183 GeV) 



Rip'l) (189 GeV) 



0.0-0.1 


52.3 


± 


2.2 


± 


2.1 


55.4 


± 


1.0 


± 


1.9 


53.7 


± 


0.6 


± 


1.7 


0.1-0.2 


39.8 


± 


1.6 


± 


1.7 


42.4 


± 


0.8 


± 


1.3 


43.4 


± 


0.5 


± 


1.2 


0.2-0.3 


31.2 


± 


1.3 


± 


1.1 


32.1 


± 


0.7 


± 


0.9 


32.7 


± 


0.4 


± 


0.7 


0.3-0.4 


23.9 


± 


1.1 


± 


1.0 


24.2 


± 


0.6 


± 


0.6 


24.1 


± 


0.3 


± 


0.4 


0.4-0.5 


18.6 


± 


1.1 


± 


0.7 


19.1 


± 


0.5 


± 


0.5 


17.9 


± 


0.3 


± 


0.3 


0.5-0.6 


16.3 


± 


1.0 


± 


0.6 


13.7 


± 


0.4 


± 


0.4 


14.8 


± 


0.3 


± 


0.3 


0.6-0.7 


10.9 


± 


0.9 


± 


0.4 


11.4 


± 


0.4 


± 


0.3 


11.1 


± 


0.2 


± 


0.2 


0.7-0.8 


8.43 


± 


0.70 


± 


0.37 


9.40 


± 


0.36 


± 


0.21 


8.65 


± 


0.20 


± 


0.17 


0.8-0.9 


7.62 


± 


0.74 


± 


0.40 


7.09 


± 


0.31 


± 


0.18 


7.46 


± 


0.19 


± 


0.13 


0.9-1.0 


5.67 


± 


0.60 


± 


0.38 


5.92 


± 


0.30 


± 


0.22 


6.42 


± 


0.18 


± 


0.19 


1.0-1.2 


3.88 


± 


0.39 


± 


0.40 


4.52 


± 


0.19 


± 


0.23 


4.95 


± 


0.12 


± 


0.16 


1.2-1.4 


2.82 


± 


0.35 


± 


0.39 


3.95 


± 


0.19 


± 


0.25 


3.56 


± 


0.10 


± 


0.13 


1.4-1.6 


2.28 


± 


0.31 


± 


0.31 


2.65 


± 


0.15 


± 


0.17 


2.94 


± 


0.09 


± 


0.08 


1.6-2.0 


1.78 


± 


0.22 


± 


0.19 


1.89 


± 


0.10 


± 


0.12 


1.95 


± 


0.06 


± 


0.06 


2.0-2.5 


1.00 


± 


0.14 


± 


0.11 


1.10 


± 


0.07 


± 


0.06 


1.16 


± 


0.04 


± 


0.05 


2.5-3.0 


0.47 


± 


0.09 


± 


0.06 


0.75 


± 


0.06 


± 


0.05 


0.72 


± 


0.03 


± 


0.03 


3.0-3.5 


0.38 


± 


0.09 


± 


0.07 


0.49 


± 


0.04 


± 


0.05 


0.55 


± 


0.03 


± 


0.03 


3.5-4.0 


0.27 


± 


0.08 


± 


0.06 


0.29 


± 


0.04 


± 


0.03 


0.37 


± 


0.02 


± 


0.02 


4.0-5.0 


0.17 


± 


0.04 


± 


0.05 


0.22 


± 


0.02 


± 


0.03 


0.23 


± 


0.01 


± 


0.02 


5.0-6.0 


0.06 


± 


0.03 


± 


0.02 


0.13 


± 


0.02 


± 


0.02 


0.12 


± 


0.01 


± 


0.01 


6.0-7.0 


0.07 


± 


0.03 


± 


0.03 


0.05 


± 


0.01 


± 


0.02 


0.08 


± 


0.01 


± 


0.01 


7.0-8.0 


0.06 


± 


0.03 


± 


0.03 


0.02 


± 


0.01 


± 


0.01 


0.04 


± 


0.01 


± 


0.01 



(pT) 



0.593 ± 0.018 ± 0.016 



0.622 ± 0.008 ± 0.009 



0.647 ± 0.005 ± 0.011 



Table 10: Measured values for the momentum spectra R(p™) — 1/a ■ da^/dp 1 ^ in the 
event plane at y/s = 172, 183 and 189 GeV. The first error is statistical, the second 
systematic. The mean values are also shown. 
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out 

1 _L 


R{p 


out 

1 

_L 


) (172 GeV) 


Rip 


out 

1 

_l_ 


) (183 GeV) 


Rip 


out 

1 

_L 


) (189 GeV) 


0.0-0.1 


71.5 


± 


9 3 




9 "i 


77.2 


± 


1 9 




1 7 


76.2 


± 


n 7 

U. 1 


in 


1.0 


0.1-0.2 


59.1 


± 


2 n 


-|- 


1 7 

.1- • I 


59.1 


± 


1 


-|- 


1.1 


59.8 


± 


n 6 


-|- 


1.1 


0.2-0.3 


39.7 


± 


1.6 


± 


1.1 


41.1 


± 


0.7 


± 


0.6 


41.7 


± 


0.5 


± 


0.7 


0.3-0.4 


26.2 


± 


1.3 


± 


0.9 


27.9 


± 


0.6 


± 


0.5 


27.6 


± 


0.4 


± 


0.4 


0.4-0.5 


18.4 


± 


1.2 


± 


0.7 


19.0 


± 


0.5 


± 


0.5 


18.8 


± 


0.3 


± 


0.4 


0.5-0.6 


11.8 


± 


1.0 


± 


0.5 


11.6 


± 


0.4 


± 


0.3 


12.2 


± 


0.3 


± 


0.3 


0.6-0.7 


7.63 


± 


0.75 


± 


0.51 


8.29 


± 


0.36 


± 


0.29 


8.45 


± 


0.22 


± 


0.27 


0.7-0.8 


3.94 


± 


0.57 


± 


0.41 


5.30 


± 


0.29 


± 


0.19 


5.65 


± 


0.18 


± 


0.19 


0.8-0.9 


3.56 


± 


0.49 


± 


0.59 


4.15 


± 


0.27 


± 


0.25 


4.11 


± 


0.16 


± 


0.13 


0.9-1.0 


2.46 


± 


0.45 


± 


0.49 


2.82 


± 


0.22 


± 


0.20 


2.91 


± 


0.13 


± 


0.10 


1.0-1.2 


1.76 


± 


0.31 


± 


0.36 


1.96 


± 


0.15 


± 


0.16 


2.04 


± 


0.09 


± 


0.09 


1.2-1.4 


0.97 


± 


0.20 


± 


0.18 


1.25 


± 


0.12 


± 


0.19 


1.23 


± 


0.07 


± 


0.08 


1.4-1.6 


0.62 


± 


0.18 


± 


0.13 


0.62 


± 


0.09 


± 


0.13 


0.69 


± 


0.05 


± 


0.06 


1.6-2.0 


0.33 


± 


0.11 


± 


0.90 


0.40 


± 


0.05 


± 


0.13 


0.43 


± 


0.04 


± 


0.05 


2.0-2.4 


0.01 


± 


0.04 


± 


0.03 


0.16 


± 


0.04 


± 


0.06 


0.20 


± 


0.03 


± 


0.03 


2.4-2.8 


0.04 


± 


0.05 


± 


0.15 


0.04 


± 


0.03 


± 


0.03 


0.11 


± 


0.02 


± 


0.03 


2.8-3.2 












0.04 


± 


0.03 


± 


0.10 


0.07 


± 


0.02 


± 


0.04 


3.2-3.6 












0.01 


± 


0.02 


± 


0.01 


0.04 


± 


0.02 


± 


0.05 


3.6-4.0 






















0.02 


± 


0.02 


± 


0.03 


(pT) 


0.282 ± 0.007 ± 0.022 


0.289 ± 0.003 ± 0.008 


0.300 ± 0.003 ± 0.004 



Table 11: Measured values for the momentum spectra i?(p° ut ) — 1/cr • da ch /chp° ut out of 
the event plane at y^s = 172, 183 and 189 GeV. The first error is statistical, the second 
systematic. The mean values are also shown. 
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R(y) (172 GeV) 



R(y) (183 GeV) 



R(y) (189 GeV) 



0.00-0.33 


7.56 


± 


0.60 


± 


0.51 


7.40 


± 


0.29 


± 


0.35 


7.44 


± 


0.18 


± 


0.30 


U.oo-U.O 1 


D.ob 


± 


0.60 


± 


U.41 


7 A K 

1 .40 


± 


0.29 


± 


U.oz 


Q HQ 


± 


0.18 


± 


n on 

u.zy 


U.D (-1.UU 


6.99 


± 


0.57 


± 


0.44 


7.60 


± 


0.28 


± 


0.26 


8.31 


± 


0.17 


± 


0.25 


1.00-1.33 


6.98 


± 


0.49 


± 


0.35 


7.51 


± 


0.25 


± 


0.23 


7.95 


± 


0.15 


± 


0.19 


1.33-1.67 


6.47 


± 


0.47 


± 


0.30 


8.07 


± 


0.24 


± 


0.22 


7.79 


± 


0.14 


± 


0.19 


1.67-2.00 


7.74 


± 


0.43 


± 


0.34 


7.30 


± 


0.21 


± 


0.19 


7.36 


± 


0.12 


± 


0.16 


2.00-2.33 


6.72 


± 


0.41 


± 


0.28 


6.99 


± 


0.19 


± 


0.21 


7.02 


± 


0.11 


± 


0.14 


2.33-2.67 


6.40 


± 


0.36 


± 


0.27 


6.45 


± 


0.17 


± 


0.18 


6.40 


± 


0.10 


± 


0.12 


2.67-3.00 


5.85 


± 


0.34 


± 


0.24 


5.89 


± 


0.15 


± 


0.16 


5.77 


± 


0.09 


± 


0.13 


3.00-3.33 


4.45 


± 


0.28 


± 


0.20 


4.96 


± 


0.14 


± 


0.13 


4.91 


± 


0.08 


± 


0.14 


3.33-3.67 


4.37 


± 


0.32 


± 


0.20 


4.00 


± 


0.13 


± 


0.13 


3.80 


± 


0.07 


± 


0.11 


3.67-4.00 


2.38 


± 


0.21 


± 


0.11 


2.97 


± 


0.11 


± 


0.12 


2.72 


± 


0.06 


± 


0.08 


4.00-4.33 


1.62 


± 


0.16 


± 


0.11 


1.76 


± 


0.09 


± 


0.08 


1.77 


± 


0.05 


± 


0.07 


4.33-4.67 


1.19 


± 


0.15 


± 


0.12 


1.00 


± 


0.06 


± 


0.07 


0.89 


± 


0.03 


± 


0.07 


4.67-5.00 


0.31 


± 


0.08 


± 


0.05 


0.51 


± 


0.05 


± 


0.05 


0.50 


± 


0.03 


± 


0.05 


5.00-5.33 


0.19 


± 


0.06 


± 


0.06 


0.22 


± 


0.03 


± 


0.03 


0.21 


± 


0.02 


± 


0.02 


5.33-5.67 


0.045 


± 


0.026 ± 


0.028 


0.089 


± 


0.017 ± 


0.015 


0.082 


± 


0.010 ± 


0.010 


5.67-6.00 


0.019 


± 


0.019 ± 


0.017 


0.039 


± 


0.011 ± 


0.009 


0.031 


± 


0.006 ± 


0.006 


6.00-6.33 


0.012 


± 


0.012 ± 


0.013 


0.007 


± 


0.004 ± 


0.001 


0.007 


± 


0.003 ± 


0.002 


(y) 


1.877 ± 0.024 ± 0.032 


1.877 ± 0.011 ± 0.026 


1.837 ± 0.006 ± 0.018 



Table 12: Measured values for the rapidity R(y) = 1/er • da^/dy distributions at y/s = 
172, 183 and 189 GeV. The first error is statistical, the second systematic. The mean 
values are also shown. 
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R(x p ) (172 GeV) 


R(xp 


) (183 GeV) 


R(x p ) (189 GeV) 


n nn n m 

U.UU-U.Ul 


r-j a r\ i on 

749. ± 30. 


± 


30. 


819 4- 
OlZ. It 


15. 


± 


21. 


855. ± 9. 


± 


20. 


U.U1-U.UZ 


A TO 1 OA 

472. ± 20. 


1 

± 


22. 


^1 K 4_ 
010. ± 


10. 


1 

± 


13. 


506. ± 6. 


1 

± 


10. 


n no n n*3 


281. ± 15. 


± 


11. 


9Q^ 4- 

zyo. it 


7. 


± 


7. 


295. ± 4. 


± 


7. 


n n 1 ? n n/i 

U.UO-U.U4 


194. ± 12. 




8. 


1 Q1 4- 

iyi. ± 


5. 




6. 


192. ± 3. 


± 


4. 


U.U4-U.U0 


loll. ± 1U. 


i 

± 


0. 


1/17 4- 
14 / . it 


A 

4. 


i 

± 


4. 


i /in 4- q 
14U. ± 6. 


i 

± 


Q 

o. 


D D^-D D6 


117. ± 8. 


± 


7. 


1 rn + 

± \Jkj. _l_ 


4. 


± 


3. 


106. ± 2. 


± 


2. 


06-0 07 


77.4 ± 7.3 


± 


4.3 


87.9 ± 


3.4 


± 


2.3 


84.4 ± 2.0 


± 


2.6 


0.07-0.08 


66.8 ± 6.0 


± 


4.0 


64.3 ± 


2.9 


± 


1.6 


66.2 ± 1.7 


± 


2.0 


0.08-0.09 


50.8 ± 5.4 


± 


3.4 


55.1 ± 


2.6 


± 


1.6 


56.0 ± 1.6 


± 


1.9 


0.09-0.10 


42.9 ± 4.9 


± 


3.2 


46.2 ± 


2.3 


± 


1.4 


46.4 ± 1.4 


± 


1.2 


0.10-0.12 


38.4 ± 3.6 


± 


2.8 


35.2 ± 


1.5 


± 


1.3 


37.0 ± 0.9 


± 


1.3 


0.12-0.14 


26.5 ± 3.1 


± 


1.5 


30.1 ± 


1.4 


± 


1.0 


27.3 ± 0.8 


± 


1.4 


0.14-0.16 


25.9 ± 2.5 


± 


1.5 


22.7 ± 


1.2 


± 


1.1 


21.1 ± 0.7 


± 


1.1 


0.16-0.18 


20.0 ± 2.5 


± 


1.3 


18.9 ± 


1.1 


± 


1.3 


16.0 ± 0.6 


± 


0.9 


0.18-0.20 


14.7 ± 2.0 


± 


1.1 


11.7 ± 


0.9 


± 


0.9 


12.1 ± 0.5 


± 


0.6 


0.20-0.25 


8.54 ± 0.93 


± 


1.11 


9.66 ± 


0.50 


± 


0.67 


9.03 ± 0.28 


± 


0.38 


0.25-0.30 


3.92 ± 0.70 


± 


0.55 


5.53 ± 


0.37 


± 


0.38 


5.20 ± 0.21 


± 


0.22 


0.30-0.40 


2.69 ± 0.34 


± 


0.37 


2.65 ± 


0.18 


± 


0.16 


2.60 ± 0.10 


± 


0.10 


0.40-0.50 


1.09 ± 0.24 


± 


0.14 


1.34 ± 


0.12 


± 


0.14 


1.13 ± 0.07 


± 


0.06 


0.50-0.60 


0.77 ± 0.18 


± 


0.12 


0.51 ± 


0.08 


± 


0.08 


0.50 ± 0.04 


± 


0.04 


0.60-0.80 


0.16 ± 0.06 


± 


0.06 


0.12 ± 


0.02 


± 


0.04 


0.14 ± 0.02 


± 


0.04 




0.0495 ± 0.0011 ± 0.0010 


0.0480 ± 0.0005 ± 0.0007 


0.0464 ± 0.0003 ± 0.0007 



Table 13: Measured values for the fragmentation functions R(x p ) = 1/er • d<7 ch /dx p at 
= 172, 183 and 189 GeV. The first error is statistical, the second systematic. The 
mean values are also shown. 
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fl(fp) (172 Gc 


V) 






(183 GeV) 




(189 GeV) 


0.0-0.2 


0.005 


± 


0.006 ± 


0.003 


0.015 


± 


0.005 ± 


0.005 


0.015 


± 


0.003 ± 


0.005 


0.2-0.4 


0.05 


± 


0.03 


± 


0.02 


0.07 


± 


0.02 


± 


0.02 


0.08 


± 


0.01 


± 


0.02 


0.4-0.6 


0.30 


± 


0.08 


± 


0.06 


0.11 


± 


0.03 


± 


0.02 


0.16 


± 


0.02 


± 


0.02 


0.6-0.8 


0.42 


± 


0.10 


± 


0.06 


0.46 


± 


0.05 


± 


0.09 


0.36 


± 


0.03 


± 


0.04 


0.8-1.0 


0.57 


± 


0.12 


± 


0.07 


0.65 


± 


0.06 


± 


0.06 


0.63 


± 


0.03 


± 


0.05 


1.0-1.2 


1.01 


± 


0.15 


± 


0.14 


1.00 


± 


0.08 


± 


0.08 


0.96 


± 


0.04 


± 


0.05 


1.2-1.4 


1.12 


± 


0.19 


± 


0.14 


1.54 


± 


0.09 


± 


0.10 


1.42 


± 


0.05 


± 


0.06 


1.4-1.6 


1.90 


± 


0.22 


± 


0.24 


2.07 


± 


0.12 


± 


0.15 


1.98 


± 


0.06 


± 


0.08 


1.6-1.8 


2.86 


± 


0.31 


± 


0.20 


2.74 


± 


0.13 


± 


0.15 


2.46 


± 


0.07 


± 


0.12 


1.8-2.0 


3.90 


± 


0.31 


± 


0.19 


3.36 


± 


0.14 


± 


0.16 


3.18 


± 


0.08 


± 


0.13 


2.0-2.2 


3.62 


± 


0.33 


± 


0.20 


3.78 


± 


0.16 


± 


0.13 


3.61 


± 


0.09 


± 


0.15 


2.2-2.4 


4.18 


± 


0.36 


± 


0.24 


4.19 


± 


0.16 


± 


0.14 


4.33 


± 


0.10 


± 


0.14 


2.4-2.6 


4.62 


± 


0.37 


± 


0.28 


4.90 


± 


0.18 


± 


0.16 


4.82 


± 


0.10 


± 


0.13 


2.6-2.8 


5.06 


± 


0.42 


± 


0.26 


5.27 


± 


0.19 


± 


0.14 


5.40 


± 


0.11 


± 


0.13 


2.8-3.0 


6.17 


± 


0.41 


± 


0.32 


5.66 


± 


0.19 


± 


0.19 


5.76 


± 


0.11 


± 


0.11 


3.0-3.2 


5.98 


± 


0.47 


± 


0.24 


6.59 


± 


0.20 


± 


0.19 


6.25 


± 


0.12 


± 


0.15 


3.2-3.4 


5.93 


± 


0.43 


± 


0.23 


6.42 


± 


0.21 


± 


0.23 


6.66 


± 


0.12 


± 


0.16 


3.4-3.6 


8.06 


± 


0.52 


± 


0.36 


6.93 


± 


0.22 


± 


0.22 


6.89 


± 


0.13 


± 


0.17 


3.6-3.8 


6.88 


± 


0.48 


± 


0.44 


7.46 


± 


0.23 


± 


0.23 


7.31 


± 


0.13 


± 


0.14 


3.8-4.0 


6.30 


± 


0.48 


± 


0.42 


7.67 


± 


0.24 


± 


0.21 


7.41 


± 


0.14 


± 


0.18 


4.0-4.2 


7.30 


± 


0.50 


± 


0.51 


7.51 


± 


0.22 


± 


0.24 


7.13 


± 


0.13 


± 


0.19 


4.2-4.4 


6.33 


± 


0.44 


± 


0.37 


7.61 


± 


0.24 


± 


0.27 


7.45 


± 


0.14 


± 


0.22 


4.4-4.6 


6.93 


± 


0.45 


± 


0.41 


6.99 


± 


0.22 


± 


0.22 


7.28 


± 


0.14 


± 


0.17 


4.6-4.8 


6.56 


± 


0.46 


± 


0.45 


6.71 


± 


0.21 


± 


0.17 


6.90 


± 


0.13 


± 


0.18 


4.8-5.0 


6.39 


± 


0.48 


± 


0.44 


6.25 


± 


0.20 


± 


0.18 


6.48 


± 


0.13 


± 


0.20 


5.0-5.2 


5.18 


± 


0.41 


± 


0.42 


5.92 


± 


0.20 


± 


0.24 


5.89 


± 


0.12 


± 


0.20 


5.2-5.4 


4.40 


± 


0.38 


± 


0.37 


4.95 


± 


0.19 


± 


0.21 


5.32 


± 


0.11 


± 


0.19 


5.4-5.6 


3.95 


± 


0.38 


± 


0.40 


4.11 


± 


0.17 


± 


0.22 


4.49 


± 


0.10 


± 


0.17 


5.6-5.8 


2.85 


± 


0.29 


± 


0.24 


3.63 


± 


0.16 


± 


0.20 


3.74 


± 


0.09 


± 


0.20 


5.8-6.0 


2.58 


± 


0.30 


± 


0.29 


2.59 


± 


0.14 


± 


0.20 


3.02 


± 


0.09 


± 


0.17 


6.0-6.2 


2.12 


± 


0.29 


± 


0.29 


2.17 


± 


0.13 


± 


0.17 


2.34 


± 


0.08 


± 


0.17 


6.2-6.4 


0.92 


± 


0.29 


± 


0.43 


1.47 


± 


0.14 


± 


0.28 


1.62 


± 


0.08 


± 


0.19 



Table 14: Measured values for the R(£ p ) = 1/cr ■ d<T c h/cl£ p distributions at y/s = 172, 
and 189 GeV. The first error is statistical, the second systematic. 



34 







172 GeV 


183 GeV 


189 GeV 


standard result 


4.031 


4.087 


4.124 


statistical error 


0.033 


0.014 


0.010 


Alternative y/s* 


0.012 


<0.001 


0.010 


Background ±5 % 


-0.002 


-0.002 


0.002 


\do\ < 5 cm 


-0.006 


0.011 


<0.001 


\z \ < 10 cm 


-0.006 


-0.006 


-0.005 


N hit9 > 80 


-0.026 


-0.023 


-0.024 


| cos6> T | < 0.7 


0.012 


-0.002 


0.002 


HERWIG 


0.005 


0.013 


0.006 


Wq C D 


0.001 


-0.001 


-0.004 


fit range 


0.024 


0.007 


0.025 


Tot. syst. error 


0.041 


0.030 


0.037 



Table 15: Results with statistical and systematic uncertainties for the position £ of the 
peak of the £ p distribution at y/s = 172, 183 and 189 GeV. 





K 


172 GeV 


183 GeV 


189 GeV 


standard result 


1.143 


1.183 


1.164 


statistical error 


0.030 


0.010 


0.008 


Alternative 


0.014 


-0.001 


<0.001 


Background ±5 % 


-0.001 


0.001 


-0.001 


\d \ < 5 cm 


-0.013 


-0.002 


<0.001 


\zq\ < 10 cm 


-0.012 


-0.012 


-0.012 


N hit9 > 80 


-0.001 


<0.001 


<0.001 


| cos6> T | < 0.7 


-0.009 


-0.006 


-0.005 


HERWIG 


0.012 


0.011 


0.005 


Wqcd 


0.005 


-0.003 


-0.007 


fit range 


0.049 


0.015 


0.026 


Tot. syst. error 


0.056 


0.023 


0.030 



Table 16: Results for the normalisation factor K in the MLLA description of the £ p 
distribution with statistical and systematic uncertainties at = 172, 183 and 189 GeV. 
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£o - (Q 


172 GeV 


183 GeV 


189 GeV 


standard result 


0.054 


0.012 


-0.035 


statistical error 


0.023 


0.010 


0.007 


Alternative \/ s' 


0.005 


-0.001 


0.008 


Background ±5 % 


0.001 


0.002 


-0.003 


\d \ < 5 cm 


0.003 


0.001 


0.004 


2 < 10 cm 


-0.002 


-0.002 


-0.003 


N hits > 80 


-0.015 


-0.015 


-0.018 


| cos6» T | < 0.7 


-0.017 


0.002 


0.006 


HERWIG 


-0.046 


-0.059 


-0.069 


Wqcd 


0.016 


0.005 


-0.009 


fit range 


0.032 


0.020 


0.020 


Tot. syst. error 


0.063 


0.064 


0.075 



Table 17: Results with statistical and systematic uncertainties for the difference of the 
position £ of the peak of the £ p distribution and its mean value (£ p )at = 172, 183 and 
189 GeV. 







172 GeV 


183 GeV 


189 GeV 


standard result 


-0.38 


-0.06 


0.16 


statistical error 


0.56 


0.44 


0.33 


Alternative 


-0.05 


0.01 


-0.03 


Background ±5 % 


<0.01 


-0.01 


0.01 


\do\ < 5 cm 


-0.03 


-0.01 


-0.02 


2 < 10 cm 


0.02 


0.01 


0.01 


N hits >S0 


0.12 


0.08 


0.07 


| cos6> T | < 0.7 


0.15 


-0.01 


-0.02 


HERWIG 


0.34 


0.26 


0.20 


Wqcd 


-0.14 


-0.03 


0.04 


fit range 


0.03 


0.02 


0.02 


Tot. syst. error 


0.42 


0.27 


0.22 



Table 18: Results with statistical and systematic uncertainties for the ratio of the 
difference of the position £ of the peak of the £ p distribution and its mean value (£ p )and the 
difference between the median and the mean of the £ p distribution, (£ — (£p))/(£m — (£p)) 
at = 172, 183 and 189 GeV. 



36 



09 10- 



1 1 1 1 1 1 1 1 1 1 p 

(a) • OPAL (189 GeV) 



~i 1 1 1 1 r 




rejected <r 



All 

Non-radiative (Z/y)* 
4-fermion backgr V///A 



_i i i i i i_ 




80 



100 





250 




200 






= 


150 


> 






100 




50 








— i — i — i — i — i — i — i — i — i — i — i — r 

rejected <; 



120 140 160 180 200 

Vs' (GeV) 

i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r 

• OPAL (189 GeV) 

(b) 

) 

T 





QCD 



Figure 1: a) Distribution of reconstructed y/s' for the data (full points) with statistical 
errors. The PYTHIA predictions for all (Z /^)* events (solid line) and for the non- 
radiative events, ^fs — \J s[ ruc < 1 GeV, (dotted line) are also shown. The hatched area 
indicates the 4-fermion background predicted by the GRC4F Monte Carlo, 
b) Distribution of the QCD event weight Wqcd- The notation is identical to that of a). 
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Figure 2: The effect of the Wqcd selection cut on the thrust distribution (see Figure |l] 
for notation). In (a) the thrust distribution, on detector level, is shown before the Wqcd 
selection is applied. In (b) the same, after the selection on M^qcd is applied. 
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Figure 3: Distributions at y/s = 189 GeV of the event shape observables thrust T, thrust 
major T ma j or , thrust minor T minor , aplanarity A, C-parameter C and heavy jet mass M H . 
Experimental statistical errors are delimited by the inner small horizontal bars. The 
total errors are shown by the outer horizontal error bars. Hadron level predictions from 
PYTHIA, HERWIG and ARIADNE are also shown. 
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Figure 4: Distributions at = 189 GeV of the event shape observables sphericity 
S, oblateness O, total jet broadening B T , wide jet broadening B w , and the transition 
value between 2- and 3-jets y^. Experimental statistical errors are delimited by the inner 
small horizontal bars. The total errors are shown by the outer error bars. Hadron level 
predictions from PYTHIA, HERWIG and ARIADNE are also shown. 
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Figure 5: Distribution of the mean thrust (T) as a function of the centre-of-mass energy, 
compared to several QCD models. 



41 




Figure 6: Distributions of the event shape observables thrust 1—T, heavy jet mass M#, 
total B T and wide B w jet broadening, C-parameter and the transition value between 2- 
and 3-jets y£, are shown together with fits of the (9(o;^)+NLLA QCD predictions, with 
a s as fitted parameter. The fitted regions are indicated by the arrows. 
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Figure 7: Fitted values of a s for each event shape at 172, 183 and 189 GeV. The inner 
error bar is statistical, the outer corresponds to the total systematic uncertainty. The 
light band correspons to the statistical uncertainty of a s obtained as the weighted mean. 
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Figure 8: Values of a s , as determined from fits to 1 — T, Mh, Bt and B\y, as function 
of energy, on a linear scale. The curves show the 0(a^) QCD prediction for a s (Q) using 
a s {Mip) = 0.119 ± 0.004 [f5|; the full line shows the central value while the dotted lines 
indicate the variation given by the uncertainty. 
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Figure 9: Values of a s as function of energy p6| . The labels NLO and NNLO refer to the 
order of calculation used. NLO corresponds to 0(a 2 s ) in e + e~ annihilations, and NNLO 
to 0(o? s ). The label Lattice refers to a s values determined from lattice QCD calculations. 
The curves show the 0(a 3 s ) QCD prediction for a s (Q) using a s (M z o) = 0.119 ± 0.004; 
the full line shows the central value while the dotted lines indicate the variation given by 
the uncertainty. 
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Figure 10: (a) Corrected distribution of the charged particle multiplicity n c h- (b) Mean 
charged multiplicity measurements by LEP experiments over a range of \/s from 91.2 GeV 
to 189 GeV. The NLLA QCD prediction for the evolution of charged particle multiplicity 
with y/s uses all data points available from 12 GeV up to 161 GeV. Also shown are the 
predictions from PYTHIA, HERWIG, and COJETS. In (b) the curve for the ARIADNE 
prediction is almost indistinguishable from the PYTHIA prediction and is omitted. Note 
that the systematic uncertainties on the OPAL measurements at energies \fs > 91.2 GeV 
are highly correlated. 
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Figure 11: (a,b,c) Distributions of the momenta p™, p° ut and of the rapidity, y; (d) shows 
the fragmentation function 1/a ■ da^/dxp with x p = 2p/ ' \fs, all compared to PYTHIA, 
HERWIG, COJETS and ARIADNE predictions. 
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Figure 12: (a) Distribution of £ p = ln(l/x p ) for charged particles at y/s— 189 GeV. 
Also shown are a fit of a skewed Gaussian and predictions by MLLA QCD, PYTHIA and 
COJETS. The curve for the ARIADNE prediction is almost indistinguishable from the 
PYTHIA prediction and is omitted, (b) Evolution of the position of the peak of the £ p 
distribution, £ , with cm. energy y/s, compared with a fit of a MLLA QCD prediction to 
the previous measurements and with predictions by PYTHIA, HERWIG and COJETS. 
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Figure 13: Results for the normalisation factor K in the MLLA description of the £ p 
distribution, compared with results from fits to the £ p distribution at lower energies |54| . 
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